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Abstract. Understanding the observer-dependent nature of quantum entanglement 
I has been a central question in relativistic quantum information. In this paper we 

will review key results on relativistic entanglement in flat and curved spacctimc and 
" pH . discuss recent work which shows that motion and gravity have observable effects on 

I entanglement between localized systems. 

1. Introduction 

m 

^ I In quantum information non-classical properties such as entanglement are exploited to 

I improve information tasks. A prototypical example of this is quantum teleportation 

O I where two observers Alice and Bob use two quantum systems in an entangled state 

^ I to transmit information about the state of a third system. Impressively, cutting-edge 

experiments involving entanglement based communications are reaching regimes where 
relativistic effects can no longer be neglected. Such is the case of protocols which 
^1 involve distributing entanglement over hundreds of kilometers [1, 2]. Understanding 

\ entanglement in relativistic settings has been a key question in relativistic quantum 

information. Early results show that entanglement is observer-dependent [3,4,5,6]. The 
entanglement between two field modes is degraded by the Unruh effect when observers 
are in uniform acceleration. We also learned that the spatial degrees of freedom of global 
fields are entangled, including the vacuum state [7, 8, 9, 10]. This entanglement can 
be extracted by point-like systems and in principal be used for quantum information 
processing (see for example [11, 12, 13, 14]). Most of the early studies on relativistic 
entanglement in non-inertial frames involved global modes. However, more recently, 
researchers in the field have focused their attention on understanding entanglement 
between fields or systems which are localized in space and time. The motivation for 
this is that entangled localized systems can be in principle measured, transformed and 
exploited for quantum information tasks. Among the most popular systems considered 
for this purpose are moving cavities [15, 16, 17, 18, 19, 20, 21], point-like detectors 
[11, 22, 23, 24, 25] and localized wave-packets [26, 27, 28]. In this paper we will review 
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global mode entanglement in flat and curved spacetime which constitutes the first step 
in the study of entanglement in quantum field theory. We will then discuss more recent 
ideas on entanglement which show that motion and gravity have observable effects on 
quantum correlations between localized systems [16, 17, 18, 19, 20, 21]. Interestingly, 
in these settings it is possible to generate quantum gates through motion in spacetime 
[19, 21, 29]. 

The observer-dependent nature of entanglement is a consequence of the particle 
content being different for different observers in quantum field theory [30, 31]. In 
fiat spacetime, all inertial observers agree on particle number and therefore, on 
entanglement. Entanglement is well defined in that case since inertial observers play 
a special role. However, in the case of curved spacetime, the entanglement in a given 
state varies even for inertial observers (see discussion in [31]). 

In special relativity one also finds that quantum correlations are observer- 
dependent. The entanglement between two spin particles is invariant only when 
the spin and momentum of the particles are considered to be a single subsystem. 
If only spin degrees of freedom are considered, different inertial observers would 
disagree on the amount of the entanglement between the particles. Some works show 
that spin entanglement in transformed into momentum entanglement under Lorentz 
transformations while some recent papers argue that considering spin degrees of freedom 
alone (by tracing over momentum) lead to inconsistencies (this will be discussed further 
in section 6). 

The paper is organized as follows: in the section (2) we will introduce technical 
tools in quantum field theory and quantum information. We will review the basics 
of field quantization focusing on the free bosonic massless case. We will describe 
the interaction of the field with point-like systems better known as Unruh-DeWitt 
detectors. By imposing boundary conditions we will describe fields contained within 
moving mirrors (cavities) and show how to construct wave-packets that are localized in 
space and time. A brief discussion on fields in curved spacetimes will be presented. We 
will end the section by reviewing measures of entanglement in the pure and mixed case 
as well as introduce the covariant matrix formalisms which allows for relatively simple 
entanglement computations in quantum field theory. In section (3) of this paper we will 
review the results on free mode entanglement in non-inertial frames, in an expanding 
universe and in a black hole spacetime. We will present ideas on how to extract field 
entanglement using Unruh-DeWitt detectors in section (4). In section (5) we will present 
a more modern view on the study of entanglement in quantum field theory where the 
entanglement between the modes of moving cavities is analyzed and review recent 
work on how localized wave-packets can be used to implement quantum information 
protocols. For completeness, in section (6) we review the concept of observer dependent 
entanglement for the case of zero acceleration. Here we discuss the Wigner rotation, 
the change in state under Lorentz transformations and their effect on entanglement for 
spin i particles and photons. Finally, in section (7) we will point out open questions, 
discuss work in progress and future directions in the understanding of entanglement in 
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quantum field theory. 

2. Technical tools 

2.1. Quantum field theory 

The theoretical framework in which questions of relativistic entanglement are analyzed 
is quantum field theory in fiat and curved spacetime. In the absence of a consistent 
quantum theory of gravity, quantum field theory allows the exploration of some aspects 
of the overlap of relativity and quantum theory by considering quantum fields on a 
classical spacetime. The most important lesson we have learned from quantum field 
theory is that fields are fundamental, while particles are derived notions (if at all 
possible) [30]. Field quantization is inequivalent for different observers and therefore, 
the particle content of the field may vary for different observers. For example, the 
Minkowski vacuum seen by inertial observers in fiat spacetime corresponds to a state 
populated with a thermal distribution of particles for observers in uniform acceleration 
[32]. The temperature, known as the Unruh temperature, is a function of the observer's 
acceleration. As we will see, a consequence of this is that the entanglement of free field 
modes in flat spacetime is observer- dependent [3, 4], and effects quantum information 
processing tasks such a teleporation [5, 33, 11]. Another interesting example is that of 
an expanding universe [30]. The vacuum state for observers in the asymptotic past is 
populated by particles as seen by observers in the future infinity [30] . The expansion of 
the universe creates particles and these particles are entangled [31, 34]. It might at first 
sight seam surprising that the dynamical Casimir effect is closely related to the Unruh 
effect [35, 36]. Both effects are predictions of quantum field theory. The vacuum state 
of an inertial cavity defined by inertial observers is inequivalent to the vacuum state 
of the cavity undergoing uniform acceleration as seen by observers moving along with 
the cavity (Rindler observers) [15, 16, 17, 18, 19, 20, 21]. Therefore, if a cavity is at 
rest and the field is in the vacuum state, entangled particles will be created when the 
cavity subsequently undergoes non-uniform accelerated motion [16, 18]. Related to this 
effect is the dynamical Casimir effect where the mirrors of the cavity oscillate [36, 37]. 
Before we discuss in more detail the entanglement between the modes of a quantum 
field in these and other scenarios we will revisit basic concepts in quantum field theory, 
considering the simplest case: the massless uncharged bosonic field (which we denote (j) 
) in a flat (l-l-l)-dimensional spacetime. Throughout our paper we will work in natural 
units c = h= 1 and the signature of the metric (-I-, — ). 

2.1.1. Global fields The massless real bosonic quantum field obeys the Klein-Gordon 
equation 00 = 0, where the d'Alambertian operator □ is defined as 




(1) 
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where g = det{gab) and 9^ = In flat (1 + l)-dimensional spacetime the metric 

is g^iu = Vfiv = {H — } and thus, □</> = d'f — d^. Minkowski coordinates {t,x) are a 
convenient choice for inertial observers. The solutions to the equation are plane waves 

w^,M(t,x) = ^i=e-^-(*-^^), (2) 

where the label M stands for Minkowski and e takes the value +1 for modes with positive 
momentum (right movers) and —1 for modes with negative momentum (left movers). 
The modes of frequency a; > are orthonormal with respect to the Lorentz invariant 
inner product 

(0, ^) = -t [ ird,<p - id,r)<P)d^', (3) 

where S is a spacelike hypersurface. These solutions are known as global fleld modes. 

To quantize the fleld the notion of a time-like Killing vector fleld is required. A 
Killing vector fleld K^^ is the tangent fleld to a flow induced by a transformation which 
leaves the metric invariant. This means that the Lie derivative of the metric tensor 
deflned by 

jCK9f,u = K^dxg^y + g^,\duK^ + guxd^K^, 

must vanish. When a spacetime admits such a structure it is possible to flnd a special 
basis for the solutions of = such that 

jCKUk,M = K''df,Uk,M = -iUJUk,M, 

where we have considered the action of a Lie derivative on a function. Vectors lying 
within the light cone at each point are called time-like. Therefore, if is a timelike 
Minkowski vector fleld, the Lie derivative corresponds to dt- By the action of the Lie 
derivative on the solutions of the Klein-Gordon equation we can identify the parameter 
uj > with a frequency, and classify the plane waves such that Uk,M are positive 
frequency solutions and jv/ negative frequency solutions. A few words about the 
physical signiflcance of the existence of a Killing vector fleld are in order. If a spacetime 
has as Killing vector K^, one can always flnd a coordinate system in which the metric 
is independent of one of the coordinates and the quantity E = p^K^ is constant along a 
geodesic with tangent vector p^ [38]. The quantity E can be considered as the conserved 
energy of a photon with 4-momentum p^. For static observers, i.e. those whose 4- 
velocity = dx^ /dr is proportional to the timelike Killing vector as = V{x)U'^, 
one deflnes the "redshift" factor V = (K^K^Y^"^ as the norm of the Killing vector (since 
U^U^ = 1). The frequency u of the photon measured by a static observer with 4- velocity 
is given by a; = p^U^, and hence u = E/V. A photon emitted by a static observer 
1 will be observed by a static observer 2 to have frequency U2 = CJ1V1/V2. Note that 
along the orbit of the Killing vector K^^ (not necessarily a geodesic), V is constant. For 
a general 1 + 1 spacetime with coordinates x = {x^,x^), a photon p^ = [uq, i:k{x)) of 
frequency loq > and wavevector of magnitude k{x) = loq \/—goo{x)/9i{x) (such that 
g/iuP^p'^ = 0) will be measured to have frequency ujk{x) = ooo \/ 9oo{x) \/(l ± Q!)/(l =F a) 
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with x-dependent Doppler factor a = \/—goo (x) / gi{x) {K^ (x) / K'^ (x) ) by a static 
observer along the orbit of the Kilhng vector K = K^{x)dxO + K^{x)dxi- In particular, 
in flat Minkowski spacetime with metric g^u = (+, — ) in (t,x) coordinates a photon of 
frequency Uq emitted by an inertial Minkowski observer will be measured to have the 
frequency ujk{x) = 00 ± a)/(l =F a) with Doppler factor a = K^{x)/K^{x). 

If the metric is static {dog^y = and goi, = 0) then the metric components are 
independent of the time coordinates t and the Klein-Gordon equation can be separated 
into space and time components as fu)it,x) = e~'"^^f^{x) (here {t,x) are general 1 + 1 
coordinates). The modes {f^^, f*) form a basis of the wave equation from which to define 
the notion of particles. By definition, a detector measures the proper time r along its 
trajectory. If the detectors's trajectory follows the orbit of the Killing field (i.e. the 
static observers defined above) the proper time will be proportional to the Killing time 
t. Modes that are positive frequency with respect to this Killing vector serve as a 
natural basis for describing the Fock space of particles [38]. Most importantly, under 
Lorentz transformations, timelike vectors are transformed into timelike vectors, so that 
the separation of modes into positive and negative frequencies remains invariant under 
boosts. In a general curved spacetime, the non-existence of a Killing field implies that 
the separation of modes into positive and negative frequencies is different along each 
point of the detectors's trajectory, and hence the concept of "particle" is lost (for further 
details, see [30] and Chap. 9 of [38]). Note that the photons of measured frequency 
ujk{x) in the previous paragraph are not pure plane waves along the Killing orbit, and 
therefore must be decomposed into the natural positive and negative frequency modes 

( ful ! fuj ) • 

Having identified positive and negative modes, the quantized field satisfying Dcp = 
is then given by the following operator value function 

Uk,Mak,M + u\]^ja\j^j)dk, 

where the creation and annihilation Minkowski operators a\ and ak^M satisfy the 
commutation relations [(^i M^^k',M] = ^k,k'- Note that the solutions have been treated 
differently by associating creation and annihilation operators with negative and positive 
frequency modes, respectively. The vacuum state is defined by the equation aA:,Af|0)"^ = 
and can be written as lO)"^ = |Ofc)"^ where 10^)"^ is the ground state of mode k. 
Particle states are constructed by the action of creation operators on the vacuum state 

K,...,n,)^ = (ni!,...,n,!)-V2«W...4^^|0)-M. 

Only when there exists a time-like Killing vector field it is meaningful to define particles. 
Observers flowing along timelike Killing vector fields are those who can properly describe 
particle states. This has important consequences to relativistic quantum information 
since the notion of particles (and therefore, subsystems) are indispensable to store 
information and thus, to define entanglement. However, in the most general case, curved 
spacetimes do not admit time-like Killing vector fields. 
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Interestingly, in the case where the spacetime admits a global timelike Killing vector 
field, the vector field is not necessarily unique. Consider two time-like Killing vector 
fields dr and df. It is then possible to find in each case a basis for the solutions to 
the Klein-Gordon equation {uk,ul} and {uk,ul} such that classification into positive 
and frequency solutions is possible with respect to Ot and df respectively. The field is 
equivalently quantized in both bases, therefore 



Using the inner product, one obtains a transformation between the mode solutions and 
correspondingly, between the creation and annihilation operators. 



where akk' = {uk,Uk') and f3kk' = ~{uk,u*^i) are called Bogoliubov coefficients. Since the 
vacua are given by ak\^) = afc|0) = it is possible to find a transformation between the 
states in the two bases. We note that as long as one of the Bogoliubov coefficients Pkk' 
is non-zero, and the un-barred state is the vacuum state, the state in the bared basis 
is populated with particles. Therefore, different Killing observers observe a different 
particle content in the field, i.e. particles are observer-dependent notions. 

In flat spacetime there are two kinds of observers who can meaningfully describe 
particles for all times: inertial observers who follow straight lines and observers in 
uniform acceleration who's trajectories are given by hyperbolas parameterized for 
example by 



where a is the proper acceleration at the reference worldline x = V*^ with proper 
time rj. (The notion of defining particles in a general curved spacetime is addressed in 
e.g. [30, 38]. For the other special cases when the acceleration (i) is asymptotically 
uniform in the past/future see e.g. [30, 38, 39], or (ii) asymptotically zero in the 
past but asymptotically uniform in the future and see e.g. [40]). The transformation 
suggests that a suitable choice of coordinates for uniformly accelerated observers are 
(?7, x) which are known as Rindler coordinates. The transformation in Eq. (4) is defined 
in the region |x| > t known as the (right) Rindler wedge I. When 77 — oo then 
t/x = tanh(a?7) — > 1 ^ x = t. Uniformly accelerated observers asymptotically approach 
the speed of light and are constrained to move in wedge I. Since the transformation 
does not cover all of Minkowski spacetime, one must define a second region called (left) 
Rindler wedge II by considering a coordinate transformation which differs from Eq. (4) 
by an overall sign in both coordinates. Rindler regions I and II are causally disconnected, 
and the lines x ± t = at 45 degrees define the Rindler horizon, Fig.(l). 

The metric in Rindler coordinates takes the form ds^ = {a^x^ dv'^ ^^x^) where a^x^ 
acts as an effective gravitational potential g-qriix) fo^^ the Rindler observer's local redshift 
factor. The Klein-Gordon equation in Rindler coordinates is {a~'^d'i — df, ^)(j) = and 





k' 



X = X cosh (arj) , 



t = X sinh {ar]) , 



(4) 
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Figure 1. Rindlcr space-time diagram: lines of constant position x ^-re hyperbolae 
and all curves of constant 77 are straight lines that come from the origin. An uniformly 
accelerated observer Rob travels along a hyperbola constrained to either region / or 
region //. 



the solutions [41, 42, 43] are again plane waves, though now with logarithmic spatial 
dependence Inx, (compare with (2)) 

u^,i = e^«W'^)i'^x-.) = 1 (^y'' - un,j, (5a) 

y/iiTu V47rfi \ In J 

y/Anuj V47rfi \ In J 

In the above u > 0, e = 1 corresponds to modes propagating to the right along 
lines of constant x — t, and e = — 1 to modes propagating to the left along lines of 
constant x + t. In the second equality we have introduced a positive constant /q of 
dimension length, and defined the dimensionless positive constant Q = u/a. Some 
authors [42] choose to label the Rindler mode by the (positive) frequency u, while other 
authors [41, 43] label the modes by the (positive) dimensionless quantity Q. (Note that 
—00 < eu/a = eQ < 00 acts as the effective wavevector for the Unruh modes, if one 
where to push the analogy with the inertial Minkowski modes (2)). Here we follow 
derivations from [43] and throughout this work, it will be understood that a wavevector 
subscript k on Minkowski modes {uk^M, etc. . . ) takes values in the range —00 to 00, 
while for Unruh modes {ukj, Ukji etc. . . ) it takes values from to 00. 

The solutions Ukj and u*^ ^ are identified as positive and negative frequency 
solutions, respectively, with respect to the timelike Killing vector field 9^. These 
solutions have support only in the right Rindler wedge and therefore are labeled by 
the subscript I. Note that they do not constitute a complete set of solutions. The 
transformation which defines Rindler region // also gives rise to the same spacetime. 
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However, the future-directed timelike Killing vector field which in this case is given by 
= —dr,, and the solutions are 

^^,n = ^ e^i-<^/'^^H-x)^^v) = ^ f'l^) ''''' ^ (6«) 

„* - ^ -i(-eiu;/a)ln(-x)+iOv) _ ^ f et - x \ 

VAnco V47rl] \ In J 

with support in region JJ. The solutions of region I together with the solutions in region 
II form a complete set of orthonormal solutions. Therefore, we can quantize the field in 
this basis as well, 



J (unjanj + unjianji + h.c.) dVt. 



Since region / is causally disconnected from region //, the mode operators in the 
separated wedges commute [a^j, a\^, jj\ = 0, etc. The vacuum state in the Rindler 
basis is \0)j^ = |0) ® |0) where akj |0) = and akji |0) = 0. Making use of the 
inner product we find the Bogoliubov transformations. 



Ofc.M — J {uk,M,un,i)o,n,i + {uk,M,Unj)o.lij 

where, for example, {uk,M, un,i) = J {K,m dtUn,i - (dtUnj) ul j^.^)dx. Upon computing 
the inner products [41, 43] the above formula can be written as 

/•oo 

auj,M = dQ [(aff^)*(cosh(rQ)aQ,/ - sinh(rQ)aJj jj) 
Jo 

+ ((^tnYi- ^^^Hrn)alj + cosh(ro)aQ,//)], (7) 



where 



„R (, ,r\ien (, ,j\—ien 

^^n = —fK= M) > "o^n = —fK= (^') 



are the Bogoliubov coefficients for the massless case, and / is an overall constant of 
dimension length, independent of and e. 

The Minkowski creation and annihilation operators result in an infinite sum of 
Rindler operators. An alternative basis for the inertial observers known as the Unruh 
basis can significantly simplify the transformations between inertial and uniformly 
accelerated observers. The Unruh modes a^.ij, an,L are appropriately chosen linear 
combinations of right-moving and left-moving Rindler modes respectively such that they 
are analytic across both regions J and II. That is, mq^/ and u*^ n are both proportional 
to (x — etY"^ when (— l)**^^ = (^^^■^y-<^^ = g-eT^ factored out of the latter region // 
mode. The Unruh modes are given by the direct Bogoliubov transformation with region 
/ and II Rindler modes for each value of Q as 

an,R = cosh{rn)an,i - sinh(rn)a^_^^, 

^ki = - sinh(rn)af7,/ + cosh{rn)alji, (8) 
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where tanh(rQ) = e~'^^. Here a^^R annihilates a right (R) moving Unruh mode travehng 
along lines of constant x — t in both wedges / and //, while a^^L annihilates a left 
(L) moving Unruh mode traveling along lines of constant x + t, with [an,R,ciQR] = 1, 
[0'Q,L, o-li l\ ~ 1 cross commutators vanishing. In terms of mode functions, the 

Bogoliubov transformation from the Rindler to the Unruh modes is given by 

unfi = cosh.{r)unj + sinh(r)M^^^^, 

uh,L = sinh(r)un,/ + cosh{r)u*^jj, (9) 

which are analytic inx—t across both Rindler wedges / and //. Note that the sign of the 
momentum k in region II is opposite of that in region J, but coupled with utilizing the 
complex conjugate of the region II Rindler mode, renders the resulting Unruh modes 
uq^r and right-movers (see [38] Chap. 9.5 for further details). 

The most general Unruh annihilation operator of purely positive Minkowski 
frequency is a linear combination of the two R, L Unruh creation operators, 

an,u = qLan,L + qRan,R, (10) 

where qi and Qr are complex numbers with Ig^p + jg/jp = 1. The introduction of the 
Unruh modes allows us to write the Minkowski annihilation operator (7) as a linear 
combination of only Unruh annihilation operators 

POO 

a^,M= dn[{al^^yan,R + {al^n)*(^n,L]- (11) 
Jo 

Hence, both Minkowski and Unruh annihilation operators annihilate the Minkowski 
vacuum, i.e. a^^M\0)M = 0, ctfi,R|0)x = 0, and aQ^LlO)^ = 0, and therefore, the Unruh 
vacuum and the Minkowski vacuum coincide. 

Using the direct Bogoliubov transformation between the Unruh and Rindler 
annihilation operators it is straight forward to show that (see e.g. [32, 38, 41]) 

lO'^)^ = — rTrE^^^h^W K)'l^fc)''' (12) 

cosh r ^-^ 

where tanhr = e"'^'^/". The vacuum state in the Rindler basis corresponds to a two 
mode squeezed state. Since the accelerated observer is constrained to move in region 
/ one must trace over the states in (the causally disconnected) region //. The density 
matrix of the Minkowski vacuum is given by po = |Ofc) (Ofc|^ and therefore, the state in 
region / corresponds to the following reduced density matrix 

= (1 - e-2-"/") ^(e-2™/a)n 1^)^ (^|/ ^ 

n 

which corresponds to a thermal state with temperature Tu = (where fc^ is the 

Boltzmann constant) proportional to the observer's acceleration. The temperature is 
known as the Unruh temperature. This is the well known Unruh effect [32]: the vacuum 
state as seen by inertial observers is a thermal state for observers in uniform acceleration. 
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2.1.2. Unruh-DeWitt detectors In order to give a more physical interpretation to the 
Unruh effect, Unruh-DeWitt detectors were introduced [32, 41, 42, 44]. The detectors 
consist of a point-like system endowed with an internal structure which can be either 
a two level system or a harmonic oscillator [45, 46]. The detector moves in spacetime 
following a classical trajectory given by x{t) where r is the detector's proper time. The 
detector couples to the field locally via it's monopole moment. Therefore, the interaction 
Hamiltonian for a harmonic oscillator detector is given by 

Hj(t) = J dkX{T){de'^^+d^ e-'^^) (^afce'('^"("^-"*(")) + aie-*(^"(")-*("») ,(13) 

where d and d'^ are annihilation and creation operators for the internal degrees of freedom 
of the detector, Q the frequency of the detector and w = the frequency of the field 
modes. The coupling function A(r) between the field and the detector can be chosen 
such that the interaction is switched on and off adiabatically. This removes transient 
effects. 

Standard calculations [30] consider a detector that follows either an inertial 
trajectory (say x{t) = vt for a detector with zero acceleration) or a uniformly accelerated 
one {x = {t^ + a~'^Y^'^). At time t = —tq the field is in the vacuum state |0)/ and the 
detector in it's ground state \0)d. The detector is then turned on and the transition 
probability of the detector to an exited state at time tq is calculated using perturbation 
theory. The transition rate per unit of proper time of the detector to a state with n = 1 
excitations at first order is given by 

1 r°° 

where {ip) j is the final state of the field which is not observed, and hence, averaged over 
all possible outcomes. This leads to the detector response function [30, 41] 

F{U) = hm hm — / / e-^^ir-r')-s\rl-slr' ^(^^ ^r^ ^ 

40 TO too 2to 

taking A(r) = exp(— s|r|), and defining g{T,T') = G{x{t),x{t')) where G{x,x') = 
f{0\ 0(x) 0(x') |0) / is the positive frequency Wightman function (WF), and the field 0(x) 
is given by the integral over all momentum of the rightmost term in large parentheses in 
(13). In the above, the response function F[u}) is a consequence of Fermi's golden rule 
and depends only on the field and not on the structure of the detector, since the latter 
has been completely factored out [30]. In general the WF is inversely proportional to the 
squared geodesic distance a{Ax) = g^u{x — x'Y{x — x'Y and hence is singular as x — )■ x' . 
This singular behavior is typically regularized by the "ie" prescription which treats F{(jj) 
as a contour integral in the complex Ar = t — t' plane. For an inertial trajectory x = vt 
one has [30] 5'(Ar) = — [27r(Ar — ie)]~^. For > the contour is closed in the lower 
half-plane for the integral to be convergent. However, since the pole in 5'(Ar) is in the 
upper half-plane, no contributions are obtained resulting in Fiuj) = as expected. For 
the uniformly accelerated trajectory x(r) = (^^(t) + a~^)^/^ with t(r) = sinh(ar) 
one obtains 5'(Ar) = — [(47r/a) sinh(aAr/2 — zea)]~^. A contour integration [30] now 
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picks up contributions from the poles along the negative imaginary axis, leading to a 
response function proportional to Ci;/(exp(27ra;/a) — The appearance of the Planck 
factor (exp(27rci;/a) — indicates that the equilibrium reached between the accelerated 
detector and the field in the vacuum state |0)/ is identical to the case when the 
detector remained unaccelerated, but was immersed in a bath at (Unruh) temperature 
Tu = a/{27rkB)- Transition rates for Unruh-DeWitt detectors in curved spacetime have 
been considered in [47]. We will see in the next section that such accelerated detectors 
can be used to extract entanglement from the field. 

Note that since the detector is point-like it couples with the same strength to every 
frequency field mode. However this is not a very realistic situation. A more physical 
situation corresponds to a detector which has a spatial profile [41, 48, 49]. Such detectors 
couple to a distribution of field modes which depends on the specifics of the profile. 
The use of detectors with infinitesimal, though non-zero, spatial extent can be utilized 
to give a more physical interpretation of the mathematically formal "ie" prescription 
used above in the computation of the detector response function. In brief, the use 
of infinitesimal sized detectors smears out the field (f){x) along the detector trajectory, 
thus regularizing the singular behavior of the WF as x — ?■ x'. Following Takagi (see 
§3.2 of [41]), one can consider the instantaneous non- rotating rest frame of the detector 
(the Fermi-Walker (FW) frame, see §13.6 of [50]). The coordinates transformation 
t(r, () = {a~^ + () sinh(ar) and x(r, () = {a~^ + Q cosh(ar) with metric ds^ = 
(1 + aQ'^dr'^ — d^"^, describes the observer's "local laboratory" where the detector sits 
at the origin ( = of the local spatial coordinates. For a "rigid" detector, its monopole 
field A4{x{t, ()) (the term in leftmost parentheses in (13) can be written in a factorized 
form Ai{x{T,Q) = ■^i''') ^^{^) fiO- The detector-field interaction is now given by 
Hint = J dr J d(M.{x{T, ()) <j){x{T, Q) which is now of the form A(r) M(r) 0'('r), where 
'P'i.'^) = J dCfiC) (Pi^iT, C) is the smeared out quantum field. Using the mode expansion 
0(x(r)) = J dk akUk{x{T)) with plane waves Uk{x{T)) = e^^'^^^'^'^'^ allows the smeared field 
to be written as 0'(r) = j dkakUk{x{T))fk{T) + h.c. Here /^(r) = / c;(e*^MK(-,C)-^nr)] = 
J dCe^^!^'^^^''^'^'^'^/^^^^^^ are wavepackets that can be effectively replaced with fk ~ e~'^^^^'^ 
(upon substitution of the coordinate transformations to the FW frame) where e is a small 
positive quantity on the order 1/a. This latter form of fk is just the "ie" prescription, 
and shows the advantage of utilizing smeared fields/wavepackets to regularize divergent 
quantities, as well as lending physical interpretation to mathematical procedures. Work 
in progress shows that a uniformly accelerated detector with a Gaussian spatial profile 
naturally couples to a peak distribution of Rindler modes [51]. Such a detector also 
naturally couples to distributions of Unruh modes which in this case yields a frequency 
distribution with two peaks corresponding to left and right moving Unruh modes. This 
detector model will help in the understanding, from physical perspective, the nature 
of Unruh modes which have been used to analyze the degradation of entanglement in 
non-inertial frames. 
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2.1.3. Moving cavities As mentioned previously, in order to preform quantum 
information tasks observers need to have access to the state and therefore, locahzing 
field modes in space and time becomes relevant. This can be achieved by confining 
quantum fields in cavities which can move in spacetime [15, 16, 17, 18, 19, 20, 21]. 
To confine a massless field within a cavity of length L appropriate Dirichlet boundary 
conditions must be imposed on the solutions of the Klein-Gordon equation given by (2) 
at the cavity mirrors placed at xi and Xr = xi + L. The normalized field solutions in 
this case correspond to sine functions 

Un,M = ^= sin {Ujn[x - Xi]) 6"^'^"*, (15) 

y/nir 

where n & N labels the energy states which now have a discrete spectrum given 
by Un = utt/L. The Killing vector field is dt which classifies Un,M and U* as 
positive and negative frequency solutions. The normalization of the solutions gives 
[Un,M, Um,M] = [U^My m] = ^nm and the mix products vanish. The field is therefore. 



n,M^n,M)i 

where creation and annihilation operators satisfy standard bosonic commutation 
relations. We have used capitalized letters for the cavity modes to distinguish them 
from global modes. 

We can also consider the field confined to a cavity in uniform acceleration by 
imposing uniformly accelerating boundary conditions to the modes in Eq. (5 a) at x = Xr 
and X = Xi- Here we considered Rindler coordinates (r], x) which are a suitable choice 
for this case once more. We consider the cavity to move in region / without loss of 
generality. The cavity is constructed such that it has a constant proper length as 
measured by a co-moving observer. Therefore, the mirrors move with different proper 
acceleration Ar = 1/xr and Ai = 1/xi. Taking advantage of the invariance under the 
boost Killing vector 9^ during acceleration we write 

[/„,^ = ^sin(^]„[x-X^])e-^""^ 



where Qn = nir/L' with L' = Xr — Xi = a~^ln(l -|- aL). For the case aL/c^ -C 1 we 
have L' ^ L, and the small difference of the acceleration at the ends of the mirror 
can be ignored with respect to the acceleration at the center of the cavity. Since 
we are only interested in the field within the mirrors, we dropped the index / and 
instead included R to denote the modes of the cavity, which is undergoing uniform 
acceleration. The quantum field inside the cavity as seen by a co-moving observer is 
given by 0r(?7,O = T^ni^nAv^ 00'n,R + KA^^Oo-I^r) where a^j^ and an,R are once 
more creation and annihilation operators with [an,R, al^, /j] = Sny- The ground state, in 
this case, is defined by fln.ij |0)^ = 0, Vn. Note that the Rindler coordinates completely 
cover the region inside the cavity and that the horizon always lies outside the cavity for 
all values of a. We assume the cavity's mirrors to be perfectly reflecting therefore, the 
states inside the cavity it will remain the same for all times [15]. This is in agreement 
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with Schutzhold and Unruh's [52] comment concerning how the non-transmissive cavity 
protects the state from Unruh radiation. If one prepares a state, for example a pure 
state, as long as the cavity remains either inertial or in uniform acceleration, the state 
remains unchanged, i.e. pure. 

More interesting is to consider the case in which the cavity is placed in a spaceship 
initially at rest which uniformly accelerates for a finite amount of time and finally moves 
at constant velocity, as illustrated in figure 2.1.3. This trajectory is interesting since 

ct 
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Figure 2. Non-uniformly accelerated trajectory from initial constant velocity to 
uniformly accelerated motion. 



it is more realistic (the uniformly accelerated trajectory assumes constant acceleration 
from past infinity to future infinity). Further, general travel scenarios can be crafted 
by considering sequences of it, and therefore it is called the basic building block 
trajectory [16]. The modes inside an cavity initially at rest will be affected by the 
non-uniform accelerated motion giving rise to particle creation and entanglement 
[16, 17, 18, 19, 20, 21],. It is possible to find Bogoliubov transformations between the 
modes of the inertial cavity and the modes of the cavity after any given travel scenario. 
The new modes will be given by 



0'm,F 



(16) 



where F denotes the modes in the final region, and a^n = {Um,F, Un^u) and = 
— {Um,F) U*^j^) are Bogoliubov coefficients. To study the above case in more detail it is 
convenient to work in the covariant matrix formalism which will be introduced in the 
section on quantum information. In the simple case of the building block trajectory the 
coefficients are given by 



a 



B 



B 
mn 
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L 
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F - 



ln( 
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(17a) 
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where 



mn 



mn 



L 
L 



^° dx fx 
— sin - Lo)) sin Qrn. ln(— ) 



dx sin {ujn{x — Lq)) sin I ) 



Here B stands for basic building block. These coefficients are difficult to handle both 
numerically and analytically. However, in the case where h = aL is small, the terms 
in the integrand can be expanded in a Maclaurin series such that, to second order, the 
coefficients take a simple form [16, 17, 18]. 

2.1.4- Curved spacetime For a general curved space metric g^uix) the d'Alambertian 
(see section 2.1.1) gives rises to a more complicated, and very likely non-separable, 
Klein-Gordon equation. However, in some special cases solutions can be found [30]. An 
example is a (1 + l)-dim expanding Robertson- Walker universe. This spacetime does 
not admit a global Killing vector field and therefore, it is not possible to define particles 
globally. This is also true for most spacetimes. However, the Robertson- Walker universe 
has a special property in that it is asymptotically flat in the past and future infinity 
regions where timelike Killing vector fields can be defined and employed to distinguish 
positive and negative solutions to the Klein-Gordon equation. We will consider that a 
scalar field living in this spacetime is in the vacuum state as seen by observers in the 
past infinity region, and show that the state will be populated with particles in the 
future infinity. 

The spacetime of a Robert son- Walker universe in (1 + l)-dim is given by ds'^ = 
dt^ — a'^{t)dx^ where the spatial sections of the spacetime are expanding (or contracting) 
uniformly according to the function a?{t). Considering the infinitesimal coordinate 
transformation drj = dt/a{t) and defining a^(t) = c^{t]) we obtain the metric ds"^ = 
(?{ri){drj^ — dy^). Consider c(?7) = 1 + e(l + tanh(cr?7)) where e and a are constants. 
This describes a toy model of a universe undergoing a period of smooth expansion. 
The parameter e is known as the expansion volume and a is the expansion rate. In 
the limit 77 — ?■ — oo the metric is ds^ = {drj^ — dx^) and in the limit 77 — 00 then 
ds"^ = (1 + 2e)((i?7^ — dx^)- The metric is flat in these regions where the vector field dr, 
has Killing properties. 

We now consider the massive Klein-Gordon equation (□ + m^)</) = in the 
Robertson- Walker spacetime described above where m is the mass of the field. The 
metric tensor (7^1, has components 



therefore, g = det{g^y) = —c^{ri) and the Klein-Gordon equation takes the form 




((9J - d^) + c(r/)m2)0 = 0. 



(18) 



Observer dependent entanglement 



15 



Exploiting spatial translational invariance we separate the solutions into 

so that the equation becomes 

d'nUv) + {k' + c{7])m^)U7]) = 0. (19) 

This equation can be solved in terms of two hypergeometric functions. We find two 
types of solutions [30] 

ui\v,x) = ^=^e^'™''~('^)''^^^°^^(^'')2Fi(a,/3,7i,5+), (20a) 



uf\v,x) = -=^e^'^^-^^-^'^('^)^^^^°^^(-^^Fi(a,/3,72,5-), (206) 
where the constants above are defined as 

a = l^ , /3 = , 7i = 1 , 72 = IH 5± = -(1 ±tanh(o-?7)),(21) 

a a (7 (J 2 

and the frequencies cOin = [k'^ + m^]^, Uout = [k'^ + m?{l + 2e)]^ and uj± = ^{cOout ± i^m)- 
We note that in the limit r] —oo the first solution becomes 

^/A-KUin 

while in the case rj +oo the second solution is 



ikx-ii^outV 

k /I 

One can see that in these limits the asymptotic solutions u"^^ and m^^^ to the Klein- 
Gordon equation are plane waves which can then be associated with positive mode 
solutions. The negative mode solutions correspond to u]* and n^*. Since u\ is associated 
with a plane wave at ?7 — — oo (past infinity) we call these solutions in- waves u^^^ = u^li^\ 
The solutions u'^^ = m^""*'' which are associated with plane waves at 77 — )■ +00 (future 
infinity) will be called out-waves. 

Using the linear transformation properties of hypergeometric functions one can 
write M^*"^ in terms of u^j^^^\ This is easier than calculating the Bogoliubov coefficients 
via brute force direct integration using the inner product. One obtains 

X) = aku'^r\v, x) + /^.^''^^^H^, x), (22) 



where 



^out 



Y !ii-±) r (1 — ^) ' 



\ Y (^ ii£iiL\ Y ( iiiioiit^ 

Pk= ( ^ ) \. " V'- / . (236) 

\u,J r(^)r(i + ^) ^ ' 

Here F is the Gamma functions with the properties r(l + z) = zT{z), T{1 — ix) = 
7r/[r(ix) sinh(7rx)] and |r(ix)p = 7r/[x sinh(7rx)] for z complex and x real. From 
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the above expressions we can read off the Bogohubov coefficients a^k' = Oik^kw and 
Pkk' = Pk^-kk'- Therefore the transformation between annihilation operators yields 
a™ = a^a^"* — f3la°^^\ We now consider that the state of the field in the past infinity is 
the vacuum state, 

oo 

|0)" = (g) |0)r' (24) 

fc=— oo 

and use the expression for the in-mode annihilation operator to calculate the state in 
the future infinity. Since the transformation between in and out annihilation operators 
only mixes modes of frequency k and —k, one finds that the state seen by observers in 
the future infinity is, 

|0)r = ^A^E^"l^)^»°t• (25) 

n 

where, 

sinh^(7rci;_/cr) 
sinh^(7ra;+/(7) 

The vacuum state from the perspective of observers in the remote past has particles 
in the remote future. Due to the expansion of the universe there has been particle 
creation. The state |0)*" is (once again) a two-mode squeezed state, and the particle 
creation process just described has a strong analogy with the quantum optical process 
of spontaneous parametric down conversion if one describes only the signal photons 
emerging from the end of a nonlinear crystal, while ignoring the idler photons [53, 33]. 
(Note: here the laser pump acts as as the "source" driving the expansion of the universe). 
Before reviewing work on entanglement in all of the above scenarios, we will introduce 
the basic tools to quantify it. 

2.2. Quantum information 

The main aim of quantum information is to learn how to store, process and read 
information using quantum systems. In this section we will briefly revise basic concepts 
of quantum entanglement for pure and mixed states which are considered to be key 
resources in quantum information. 

2.2.1. Entanglement Entanglement is a quantum property which is a consequence of 
the superposition principle and the tensor product structure of the Hilbert space. The 
pure bi-partite case is well understood. The state of two particles A and B is a vector 
in a {da x (i6)-dimensional Hilbert space l-iah = Ha ® T-Lb- The space l-iab is the tensor 
product of the subspaces "Ha and l-Lh of each particle. An pure state element of the space 
Hah is written as \i)ab) = Y.i,j Aj\i)a ® \ j)b- 

A state \ipab) £ Ha ®Hb is separable if \ipab) = \(p)a ® |v^)fe- Separable states can be 
prepared by local operations and classical communication. This means that observers 
manipulate each particle independently by making measurements, or applying unitary 
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transformations of the form \ipab) = Ua® Ub \(f>)a ® Iv^)^ where Ua and Ub are unitaries 
acting on particles A and B, respectively. The observers are also allowed to exchange 
classical information. If the state is not separable then it is entangled. An entangled 
state cannot be prepared by local operations and classical communication, observers 
must make global operations on the systems. To determine whether or not a general 
pure state \ipab) = J2 ^ijl"^) a\j) b is entangled we consider the following theorem: 

Schmidt decomposition: Let Ha and Tia be Hilbert spaces of dimension da 
and db, respectively. For any vector \ipab) ^ 'Ha ® Hb there exists a sets of 
orthonormal vectors {\j)a}a=i,...,da and {\l)b}b=i,...,dt '^b such that we 

can write \ipab) = J2i=i d'^iK)aK)b- The Schmidt coefficients Aj are non- 
negative scalars such that X]j -^f ~ runs over the minimum 
dimension d = min{da, db) of the two Hilbert spaces. 

This special basis is called the Schmidt basis [54, 55]. Note that the correlations 
between systems A and B are now made explicit. For example, if Aj = 1 and Xj-^i = 
then the state is separable. In the case that all the Aj's are equal and Aj = 1/ y/d then the 
state is maximally entangled. It is clear that the distribution of the Schmidt coefficients 
determine how entangled the state is. Therefore, to quantify entanglement in the pure 
bi-partite case we need a monotonous and continuous function of the Aj's such that 

(i) S'(Aj) = for separable states, 

(ii) S{Xi) = log{d) for maximally entangled states. 

Considering the density matrix pab = \'ipab){'^ab\ and it's reduced density matrix 
Pb = Tra{pab)i we find that the von Neumann entropy 

S{pb) = - Tr {pb\og2 pb) = |A,plog2 |Aip, 

i 

quantifies the entanglement between system A and B. We observe from the Schmidt 
decomposition that it is equivalent to trace over either system A or B and therefore, 
S{Pa) = S{pb)- The von Neumann entropy of a pure state is S{pab) = 0. 

Quantifying entanglement in the mixed case is more involved since now there is no 
analog to the Schmidt decomposition. However, it is still possible to define a separable 
mixed state. A bipartite mixed state is separable if we can write its density matrix as 

pab = Y^iPipi ® Pb where EiP* = 1- 

To determine whether or not a general mixed state is entangled it is convenient to 
define the partial transpose of a density matrix. Consider the general bipartite mixed 

state 

Pab = Cjjkl \i)a\j)bb{k\a{l\ = "^Cijkl \i)aa{l\ ® \j)bb{k\. (26) 

ijkl ijkl 

The partial transpose p^jf of pab is (transposing on system a) 

Pab^ = J2'^l^^^\'^^\^^>^b^^\<^(^\^ (27) 
ijkl 
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or equivalently (transposing on system b) 

ijkl 

Since the partial transpose of a separable state has positive eigenvalues it is possible 
to construct a bipartite separability criterion: If the eigenvalues of pf^^ ^ then pab 
is separable. However, this criterion is only sufficient for 2 x 2 and 2x3 systems [56]. 
For systems of higher dimension the criterion is only necessary, meaning that there are 
(bound) entangled states with positive partial transpose. Such states are known as 
bound entangled states. 

Adding up the negative eigenvalues gives an estimate of how entangled a state is. 
Therefore, we will now define the negativity and logarithmic negativity [57, 58] which 
are two entanglement monotones. The negativity of a density matrix pab is defined as 
the sum of the absolute values of the negative eigenvalues of the partial transpose p^^f, 

WPab II 



N{p. 



'ab) 



2 



where || ■ || denotes the trace norm ||X|| := Tr 



Vx^x 



a density matrix pab is defined as E]^{pab) '■= log2 \\Pab 



The logarithmic negativity of 



2.2.2. Covariant Matrix formalism The covariant matrix formalism is usually 
employed in quantum optics [59]. It is a framework that involves simple mathematical 
tools which are applicable to systems consisting of a finite number of harmonic 
oscillators. While in general computing entanglement can be very involved, it has 
been shown that in the case of Gaussian states (such as coherent, squeezed and thermal 
states) the covariant matrix formalism can be employed to produce computable measures 
of quantum information [60]. Gaussian states are described by quasi-probability 
distributions of Gaussian shape in phase space [61]. This technique is very useful in 
relativistic quantum information since Bogoliubov transformations map Gaussian states 
back onto Gaussian states [20, 62, 63]. As long as one considers initially Gaussian states 
and a finite number of field modes it is possible to calculate entanglement and other 
interesting quantities using the measures developed in this framework. In the standard 
Fock space description of quantum field theory quantum states correspond to density 
matrices which contain all the information pertaining to the state. In the covariant 
matrix formalism the notion of density matrix is replaced by a real symplectic covariant 
matrix a defined by 



cr. 



;X,X, + X,X,) - 2 (X,)(X,), (29) 



where the operators Xj are the generalized positions and momenta, i.e 

1 

3 



X(2n-1) = ^(fln + 4)) (30 a) 



X(2n) = -^(a„ - al) , (306) 
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and the index n = 1,2,3,... labels the modes, and ( O ) denotes the expectation value 
of the operator O with respect to the initial Gaussian state. The matrix a together with 
the vector of first moments ( Xj ) completely characterizes all Gaussian states. However, 
the covariance matrix is sufficient to compute entanglement. Unitary transformations 
in the Fock space are replaced in this formalism by symplectic transformations in phase 
space. A transformation S is called symplectic, if it leaves the symplectic form Q (a 



numerical matrix) invariant, i.e., S Q = Q, where = 2iQij and {Si 



The expression for the symplectic representation of a Bogoliubov transformation in 
terms of its general coefficients amn and Pmn takes a simple form [20]. The matrix S is 
decomposed into 2x2 blocks M.mn as 



S 



( A^ii 



Mx2 
M32 



Mr, 
M23 
M33 



\ 



(31) 



where the sub-blocks are given by 

n 



Mr 



-^{amn - Pmn) 5^(a^ 



+ /3. 



(32) 



Here ^{z) and '^{z) denote the real part and imaginary part of the complex number 
z respectively. The transformed covariance matrix a is then simply obtained as 
a = S a . This transformation ensures that if a is Gaussian, a remains Gaussian. 
Fortunately, for Gaussian states computable measures of entanglement have been 
developed [60, 64]. For example, the negativity between two modes is given by 

7V = max{0, (1 -P_)/2i>_}, (33) 

where z/_ is the smallest symplectic eigenvalue of the partial transpose matrix akk' = 
Tk'Ckk' Tk' , where T^/ = diag{l, 1, 1, —1} partially transposes mode k', see Ref. [61]. 

Computing the negativity between two modes k and k' involves first tracing over all 
other modes. Impressively, the partial trace over any subset of modes is computed by 
eliminating all rows and columns corresponding to all modes other than k and k'. The 
smallest symplectic eigenvalue z/_ is then obtained by diagonalizing the matrix iQ a^k' 
by a symplectic operation D, yielding the two eigenvalues < z/_ < The state akk' 
is entangled if < z/_ < 1. 



3. Entanglement of global modes 
3.1. Flat spacetime entanglement 

In the Unruh effect the Minkowski vacuum is a thermal state for uniformly accelerated 
observers. Since the Rindler regions I and II are causally disconnected, uniformly 
accelerated observers in one Rindler region (wedge) have no access to information from 
the other Rindler region. Therefore, the state of the Rindler observer, adapted to 
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his/her particular wedge by tracing out over the unaccessible wedge, is mixed. The state 
appears more mixed for observers with increasing acceleration. By this we mean that 
the temperature associated to the state by a particular observer is higher for observers 
with larger proper accelerations. This has an effect on entanglement [3, 4]. Consider a 
pure entangled Minkowski state of two field modes each of them labeled by a different 
Minkowski frequency. We assume that two inertial observers, Alice and Bob, are able to 
distinguish these modes and agree that the state is maximally entangled. As typically 
considered in the literature, Alice and a Rindler observer in region J, called Rob, also 
analyze the state. Since Rindler observers in Region I must trace over the region II part 
of the state, the correlations of the state will be partially lost. Since the vacuum itself 
exhibits greater noise for observers with higher accelerations one also expects other 
states to appear more mixed and less correlated for higher accelerations. Therefore, 
there is a degradation of entanglement due to the Unruh effect. Let us point out that 
general region / states are not necessarily thermal after tracing out over region // states, 
however, they are mixed. A thermal state is a particular type of mixed state in which 
the fraction of the ensemble in each pure state is given by a Boltzmann distribution. In 
the case of Gaussian states such as the vacuum, thermal states are obtained after tracing 
over modes. Calculations involving Minkowski states are very involved because a single 
Minkowski mode corresponds to an infinite superposition (continuous) of Rindler modes 
(7) making the trace operation highly non-trivial. It is therefore, more convenient from 
a mathematical point of view to analyze inertial states involving Unruh modes such as 
the Bell state, 

|vl>) = -L {\0k)^\0nf + \lk)^\lnf) , (34) 

where Ai and U label Minkowski and Unruh states respectively with frequencies k and 
Q = Lo/a. This state was introduced in Bruschi, et.al [43]. For more details on the 
use of Unruh modes to analyze entanglement and the single mode approximation see 
the appendix, section 8. For inertial observers the modes k and > (which must be 
two distinguishable modes) are maximally correlated since the Bell state is maximally 
entangled. It is then interesting to investigate to what degree the state is entangled 
when described by observers in uniform acceleration. In the simplest scenario shown in 
Fig. (1), Alice is considered to be inertial and a uniformly accelerated observer Rob is 
introduced, who analyzes mode VL. To study this situation, the states corresponding to 
Rob must be transformed into the appropriate basis, in this case, the Rindler basis. We 
have already calculated the transformation for the vacuum state in (8), and with that 
in hand, we find the single particle Unruh state 

4,f/|0n)" = \lnf = ^--^tanh"(r)v/^|(n + l)n),|nn),,, (35) 

cosh r ^-^ 

^ ' n 

where we have chosen \qji\ = 1 in Eq. (10). Thus, the maximally entangled state seen 
by inertial Alice and accelerated Rob is 
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+ ^^^y2tanh''{r)V^^\{n + l)n)j\nn)n. (36) 

V2 cosh(r) ^ 

Since Rob, with trajectory in region J, is causally disconnected from region // we must 
take the trace over region //. The density matrix for the Alice-Rob subsystem is 



PAR = — r27T E tanh2"(r)p„, (37) 
cosh (r) 

^ ' n=0 

where 

Pn = |Ofe,nn) (0fc,r2n| + ^ (|Ofe,nn) (U, (^ + l)n| + |lfc, {n + l)n) (Ofe,ran|) 

cosh(rj 

n + 1 

H 7^7^ (?T- + l)n) (Ifc, (?^ + l)n| , 

cosh [r) 

involves only Ai and region I states. The entanglement between the modes can be 
quantified by the logarithmic negativity yielding M{par) = log2((l/2 cosh^(r)) + S) 
where 



^ ^ tanh2"(r) // n , 2, X 4 

^-^ 2 cosh (r) y \sinh (r) / cosh (r) 

One observes that entanglement is degraded for observers in uniform acceleration 
and therefore, it is observer dependent. In the flat spacetime case, one can argue that 
this is an effect of Rob being non-inertial. Rob must be in a spaceship to be accelerated 
and energy must be supplied to accelerate the system. In the flat spacetime case, 
inertial observers play a special role and therefore, a well defined notion of entanglement 
corresponds to the entanglement described from the inertial perspective. However, in 
curved spacetime different observers describe a different particle content in the field 
which results in different degrees of entanglement in the field. In that case, there is no 
well-defined notion of entanglement. 

Using Unruh rather than Minkowski states considerably simplifies the mathematical 
techniques involved in the analysis of entanglement. However, it introduces 
interpretational issues [65]. The physical interpretation of Unruh states requires deeper 
understanding. Work in progress shows that finite size accelerated detectors naturally 
couple to peaked distributions of Unruh modes [51]. Such detectors will provide further 
insights on the states considered in Eq. (34). Currently, it is only possible to say that the 
work of Bruschi, et. al. [43] studies the entanglement of a family of maximally entangled 
states parameterized by the dimensionless Vt = u/a in non-inertial frames. By fixing 
the physical frequency uj and changing a one analyzes the entanglement in a family of 
states, all of which have the same frequency u as seen by observers with different proper 
acceleration a. An alternative view is that the entanglement corresponds to a family 
of states with different physical frequency u as seen by the same observer moving with 
fixed proper acceleration a [66]. 

The analysis of entanglement in non-inertial frames has been generalized in several 
directions and different quantum fields have been considered. For example, the work by 
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Adesso, et. al. [62] considers two-mode squeezed states in the inertial frame (instead of 
the Bell states analyzed above) allowing them to employ the covariance matrix formalism 
which is only applicable for Gaussian states. The authors studied the case when both 
observers accelerate finding that entanglement vanishes at finite acceleration. This effect 
is also known in the literature as sudden death of entanglement [67]. The effect is similar 
to that which occurs in parametric down conversion when a two-mode entangled state 
is feed into an amplification channel (see also [68]). It was shown that the rate of 
entanglement degradation depends on the entanglement in the inertial frame. States 
with higher degrees of inertial entanglement degrade more in non-inertial frames. The 
degradation of entanglement has also been analyzed for inertial thermal states and non- 
maximally entangled states [69]. 

Other types of correlations relevant to quantum information theory have also been 
studied in non-inertial frames. For example, it was found that classical correlations are 
conserved in case one observer is non-inertial [3, 4] and degraded when both observers 
accelerate [62]. More recently, the quantum discord [70, 71], a type of purely quantum 
correlations which are present even in separable states, was analyzed in non-inertial 
frames. Datta [72] found that in a regime where there is no distillable entanglement there 
is a finite amount of quantum discord. Interestingly, in the limit of infinite acceleration 
the discord does not vanish. Multi-partite entanglement has also been subject of study 
for observers in uniform acceleration [4, 62]. The bi-partite entanglement between 
modes is degraded in non-inertial frames, however, three and four-partite correlations 
are created between modes in regions / and // [62, 73]. Multipartite entanglement has 
also been analyzed when the state contains tri-partite correlations in the inertial frame 
[74]. 

The analysis of entanglement in non-inertial frames has been considered for 
different types of fields. For example, Dirac fields [4, 69, 74, 75, 76, 77, 78], the 
electromagnetic field [79] and charged bosonic fields have been studied [66]. The 
results for fermionic fields show that entanglement is degraded in non-inertial frames 
however, the entanglement remains finite in the infinite acceleration limit [4, 69, 77, 80]. 
Interestingly, for fermions the entanglement in the infinite acceleration limit does not 
violate Bell's inequalities [78]. Multipartite entanglement also presents interesting 
differences. In the bosonic case genuine tripartite entanglement is generated amongst 
the modes of Alice, Rob and Anti-Rob (an observer with uniform acceleration in region 
//) [62]. These correlations increase with acceleration. However, in the fermionic case 
tripartite entanglement is always zero [4]. Recent work [74] considers an inertial tripartite 
state which contains genuine tripartite correlations, as measured by the Svetlichny 
inequality [81]. These correlations persist for any finite acceleration for both GGHZ 
and MS three-qubit states provided the respective control parameters are appropriately 
chosen, vanishing only in the infinite acceleration limit. This indicates that tripartite 
entanglement and its associated nonlocal correlations are more robust to relativistic 
effects than their bipartite counterparts. It is important to mention that an ambiguity on 
the definition of entanglement between fermionic fields has been pointed out in [82, 83]. 
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However, these results are subject of a present debate [84, 85]. (Note that in [4], signs 
arising from the transposition of anti-commuting fermionic operators could always be 
absorbed into an explicit phase that was introduced, and hence into a redefinition of 
the fermionic creation operators). 

The main differences in the entanglement degradation for scalar and fermionic fields 
arise due to particle statistics [66, 77, 86, 87] and not to dimensionality of the Hilbert 
space. In the fermionic redistribution of entanglement between particle and 

antiparticle states seams to prevent the degradation of entanglement in the infinite 
acceleration limit [77]. Such redistribution does not occur in the bosonic case [66]. 
Initially it was conjectured that the Hilbert space dimensionality played a role [4, 69]. In 
the bosonic case the Hilbert space is infinite dimensional while in the case of Grassmann 
scalars, the lowest dimensional fermions, the fermionic Hilbert space is of dimension two. 
One might suspect that entanglement is degraded stronger in the bosonic case since for 
higher accelerations the bosonic harmonic oscillators become populated with increasing 
temperature while Grassmann scalars are always in mixed states which involve only two 
states. Considering higher dimensional fermions [88] and bosonic fields with truncated 
Hilbert spaces [86] showed that it is statistics rather than dimensionality which produces 
the observed differences in the degradation of entanglement in non-inertial frames. In 
the case of entanglement degradation for electromagnetic fields it was found that photon 
helicity entangled states do not degrade with acceleration [79]. 

3.2. Entanglement in curved spacetime. 

General curved spacetimes do not have global timelike Killing vector fields and therefore 
it is not possible to define subsystems. Without the notion of subsystem it is not possible 
to study entanglement. In spite of this, special cases have been analyzed. For instance, 
black hole spacetimes where the spacetime is approximately flat at the horizon, and 
Robertson- Walker universes which have two asymptotically fiat spacetime regions. 

3.2.1. Entanglement across black hole horizons The Schwarzchild spacetime of an 
eternal black hole describes the geometry of a spherical non-rotating mass m [30]. 
Considering only the radial component, the metric is 

ds^ = (1 - {2m/R))dT^ - (1 - {2m/R))-^dR^. 

where the horizon of the black hole is at i? = 2m. Considering the following coordinate 
change R — 2m = x'^/8m, such that 

1 - (2m R) = 8mR = , ^ , ' = , ,[ ^ iAx)^ 

when a; ~ with A = l/4m. This means that dR^ = {Ax^)dx^. Therefore, very close 
to the horizon R ^ 2m, the Schwarzschild spacetime can be approximated by Rindler 
spacetime 

ds^ = -{AxfdT^ + dx^, 
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where the acceleration parameter a = = 4m. 

This means that, very close to the horizon of the black hole, we can consider Alice 
as being an inertial observer and falling into the black hole, while Rob escapes the fall 
by being a stationary, accelerated observer. If Alice claims that the state of the field is 
the vacuum state, then Rob detects a thermal state of the form 

lO'^')^ = E l^'^)" l^'^^)"* ' (39) 

coshio) ^-^ 

where cosh(g) = (1 — e2m7rtj^-i/2 jg ^ function of the mass of the black hole m. Here in 
and out denote the modes inside and outside the black hole. \0k)^ is the state detected 
by Alice. If we then consider the case where two modes of the field in Alice's frame 
are maximally entangled, Rob will detect less entanglement between the modes due to 
the Hawking effect [3, 73, 89, 90]. Using the covariant matrix formalism an inertial 
two- mode squeezed state was analyzed [73]. The entanglement between the modes is 
degraded for observers outside the black hole, however genuine four-partite entanglement 
between the modes inside and outside the black hole are created. For low black hole 
masses entanglement vanishes and classical correlations also degrade. The analysis of 
entanglement degradation has also been considered taking into account the distance of 
the observer to the event horizon [91] by writing the metric in terms of the proper time 
Tq of a stationary observer placed at tq, 

ds^ = -{Axfdr^ + dx^ 

which yields a modified acceleration A = A/ ^/l — 2m/ro = {Amy^l — 2m/ro)~^ giving 
rise to the same state as in Eq. (39) with a modified squeezing parameter q given by 
cosh(g) = (^x_e2m7raj/-y/i-2m/ro'j-i/2_ This study has been also carried out in the fermionic 
case [87, 88], and in higher dimensional black holes [92]. 

3.2.2. Entanglement in an expanding universe Our second curved space example is 
the Robertson- Walker universe [30]. In this spacetime particles can be defined in two 
asymptotically fiat regions. Considering the field in the remote past to be in the vacuum 
state, it is found that in the distant future the state contains particles. In the interim 
region, where the universe is undergoing expansion, no sensible notion of particles exist. 
In [31] it was shown that in the future infinity region entanglement has been created 
between field modes. Interestingly, it is possible to learn about the expansion parameters 
of the Universe from the entanglement generated. The vacuum state in the past infinity 
corresponds to the following two mode states in the future infinity 

ior = v/r^5^7"in)nnr_t 

n 

Since the state is pure, we can employ the von Neumann entropy to quantify the 
entanglement generated the field modes k and —k. In order to do this we need to 
compute the reduced density matrix for one of the modes. The density matrix for the 
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state is 

po = |o)"(or" = (1 - 7) j]7^"+"'V)r»-fe(^ir*(^i-t- (40) 

n,m 

The reduced density matrix for mode k is obtained by tracing over mode —k 
pu = tr.M = (1 - l)Y.^^''^'^^Sr.m\n)r{m\r 

n,m 

= (i-7)E^'"i^)r*(^ir- (41) 

n,m 

Since the reduced denstiy matrix is aheady in diagonal form with eigenvalues A„ = 
(1 — 7)7^", it is straight forward to compute the von Neumann entropy 

r -y/l^^—i) 

S{Pk) = -(1 - 7) - 7)7'") = log2 ^ 



1-7 



(42) 



7 = ■T.r ; ^ (43) 



Entanglement has been created in the remote future due to the expansion of the universe. 
The entanglement depends on the cosmological constants since the coefficient 7 depends 
on the expansion rate a, the expansion volume e and the frequency of the modes involved 
through 

sinh^(7rci;_/o") 
sinh^(7r(X'+/(T) 

with u± = l/2{ujout ± i^in), ^in = [k"^ + ni^]^ and cOout = [k"^ + m^il + 2e)]^. In the case 
of light particles, the equations can be inverted and we can show that we can estimate 
the expansion parameters from the entanglement. For a calculation in the covariant 
matrix formalism see [20]. The analysis of entanglement in an expanding universe has 
been extended to Dirac fields [34] and entanglement in particle creation has also been 
studied in [93]. 



4. Extracting global mode entanglement using point-like detectors 

Although early investigations on relativistic entanglement allowed us to understand 
some aspects of the observer dependent nature of entanglement, the results are purely 
academic. Global mode entanglement cannot be measured nor manipulated by local 
observers. Alice would not be able to take an entangled state and apply, for example, a 
CNOT gate on it since the states live in the whole spacetime. This would require infinite 
time and energy making impossible any practical quantum information task. However, 
it has been shown that global mode entanglement can be accessed by local observers 
employing Unruh-Dewitt detectors [7, 8, 9, 10, 12]. Two Unruh-deWitt detectors at 
rest can become entangled by interacting with the Minkowski vacuum even if they are 
space-like separated. This shows that the vacuum state is entangled. The vacuum is not 
mode-wise entangled however, it contains entanglement in spatial degrees of freedom. 
This entanglement is extracted, or transferred (swapped in the quantum information 
language) from the field to the detector's degrees of freedom. In principle, it should 
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be possible to use Unruh-Dewitt detectors to study some aspects of the global mode 
entanglement analysis within a more practical approach [94]. 

Since Unruh-Dewitt detectors are local, one can consider them as suitable systems 
for the implementation of quantum information tasks. For example, Alice and Bob 
could be spacelike separated and entangle thier detectors through the interaction with 
the vacuum. This entanglement could be exploited for quantum information tasks. 
The entanglement extracted is currently too small to be useful, but in the future, 
enhancement schemes may be developed. At the moment, the main limitation in using 
Unruh-Dewitt detectors for quantum information processing is that the mathematical 
techniques involved already become un-manageable with a few systems. For example, 
for two detectors one must consider the interaction Hamiltonian 

+ J dk\{T2){d2 e^^^"^ + 4 e-*"^"2) (ak e^('=^2{r2)-<^t2{r2)) ^ ^-i(hx^{r2)-Hr2))'^ ^ 

for which entanglement calculations involve 4t/i-order perturbation theory [8, 9, 10]. 
Current research is focused on finding ways to modify the detectors in order to make 
them more mathematically accessible [51]. 

In spite of the mathematical difficulties involved, interesting results have been 
found in relativistic quantum information using point-like detectors. For example, the 
entanglement between an Unruh-DeWitt detector and a massless bosonic field has been 
considered in [22] showing that the field and detector remain entangled at late times. 
Information stored in the point-like system flows into the field propagating with the 
radiation into null infinity. The information lost by the detector is never restored (non- 
Markovian regime). Another analysis considers two detectors, one inertial belonging to 
Alice, and one in uniform acceleration belonging to Rob, which do not interact directly 
but are coupled to the same field [24]. The detectors are initially entangled and their 
entanglement vanishes after a finite time. The time at which the entanglement vanishes 
is shorter for higher accelerations, as expected. This setting has been considered recently 
for teleportation of a coherent state [11]. When both detectors are at rest and initially 
entangled, the detectors also disentangle and the dependence with their separation is 
studied in [23]. The case of Rob having a detector in non-uniform acceleration has been 
studied in [24] . The authors show that the disentangling time slows down if the motion 
is non-uniform. Two uniformly accelerated detectors have also been considered, each 
on of them in a different Rindler wedge [7, 9]. The authors show that entanglement is 
generated through their interaction with the Minkowski vacuum even if they are always 
causally disconnected. This shows again that the vacuum contains quantum correlations. 
Another interesting analysis is that of the creation and loss of entanglement between 
detectors moving in circular paths [25]. 

Recently, Olson and Ralph [10] found that by independently quantizing the past 
and future spacetime regions of a massless quantum field it is also possible to extract 
vacuum entanglement from timelike separated regions. To extract the entanglement 
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they considered two timelike separated point-like detectors at rest characterized by a 
strongly time-dependent energy gap. Time-like entanglement has been further analyzed 
in [95] and exploited in a teleportation-type protocol [13, 14] which can be implemented 
in a circuit QED experiment [14]. 

In curved spacetimes the entanglement between point-like detectors can been used 
to distinguish between two universes. Two inertial detectors in flat spacetime would 
become entangled if they interacted with a Minkowski thermal field however, they would 
not become entangled in an exponentially expanding de Sitter spacetime [96]. 

5. Accelerated cavities and localized wave-packet for quantum information 

Work on the use of Unruh-Dewitt detectors for relativistic quantum information 
processing has inspired efforts to find suitable systems to store information in the 
framework of quantum field theory, with the hope of developing mathematical techniques 
which render the calculations accessible. At the moment moving cavities and localized 
wave-packets seem plausible systems for this purpose. 

5.1. Entanglement creation and degradation in moving cavities 

The idea of using cavities moving in spacetime for relativistic quantum information 
processing was proposed in [5] where two observers Alice and Rob share an entangled 
state of two cavities as a resource for quantum teleportation in non-inertial frames. The 
mathematical techniques to implement this idea were introduced in [15] for inertial and 
uniformly accelerated cavities and in [16] for non-uniform motion. The authors in [15] 
showed that it is possible to entangle two cavities, one inertial and the second in uniform 
acceleration, by letting a single two level atom interact with the cavities. The scheme is 
a generalization of [97] in which the modes of two cavities at rest positioned next to each 
other become maximally entangled after an excited atom moving through them emits 
an excitation. Since it is not possible to know in which cavity the atom emitted, the 
cavities become entangled. When one of the cavities is in uniform acceleration the cavity 
modes become detuned affecting the atom's emission probabilities. The modes of an 
inertial and accelerated cavity are given in (15) and (16), respectively. The interaction 
Hamiltonian considered by the authors is of the Unruh-Dewitt type given in (44) only 
in the sense that the atom is a two-level ponit-like system interacting with discrete 
cavity modes. The authors show that the ability to maximally entangle the cavities 
is reduced with acceleration. However, by appropriately modifying the cavity length 
it is possible to alter the modes and compensate for the effect. An advantage of this 
scheme is that the quantum information stored in cavities is protected from the Unruh 
effect (unlike for information stored in point-like systems) by the mirrors. As long as 
the cavity walls are perfectly reflecting the fleld inside the cavity is independent from 
the fleld outside. While the cavity undergoes uniform acceleration (again we consider 
the case aL/c? <^ 1 so that the acceleration is essentially constant across the proper 
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length of the cavity), the states within the cavity however, remain pure since the Unruh 
horizon is not contained in this spacetime region [52] . The quantum information stored 
in field modes remains constant at all times as long as the cavities are inertial or in 
uniform acceleration. However the situation changes when the motion is non-uniform. 

Consider a single inertial cavity with all modes in the vacuum state. If the cavity 
accelerates for a finite time the modes will become populated as shown is section 
2.1.3. Working in the low acceleration regime /i < 1 it is possible to find analytical 
expressions for the Bogoliubov coefficients for bosons [16] and fermions [18] and calculate 
the entanglement produced between two modes k and k' using the covariant matrix 
formalism [98]. The entanglement quantified by the logarithmic negativity is given by 
A/" = \Pf^l/\ which is the first order correction to the Bogoliubov coefficient Pk^k' [20]. 
The entanglement between the modes can be calculated for any trajectory composed of 
segments of inertial and uniformly accelerated segments. The entanglement in a basic 
building block composed of initial and final inertial segments with a single period of 
uniform acceleration is periodic in the acceleration period. A consequence of this is 
that if one maximally entangles two modes each in a different cavity as proposed in 
[15] and then accelerates for a finite time, the entanglement between those modes will 
be degraded since the initial information will be transferred to other modes. Since the 
evolution of entanglement is periodic, Rob can plan his travel such that at the end of it 
his cavity will still be maximally entangled with Alice's. 

Interesting conclusions can be drawn from comparing the analysis of entanglement 
generation in fermionic and bosonic cavities [98]. In both cases entanglement is always 
generated under non-uniform motion however, the entanglement generated is different 
due to particle statistics. In the bosonic case motion can populate modes with the 
characteristic that any number of excitations can be produced in a given mode. However, 
fermionic statistics constrains the number of particles in each mode giving rise to 
differences in the entanglement produced. 

Another interesting property of motion generated entanglement is that for any 
trajectory segment constructed from periods of inertial and uniformly accelerated 
motion, there is an entanglement resonance when the segment is repeated N times [19]. 
It can be shown analytically that when the frequency associated with the segment period 
is equal to the sum of the frequencies of the two modes then the entanglement between 
those two modes grows linearly with A^. The presence of the resonance is independent 
of the segment travel details however the degree of entanglement does depend on the 
specific trajectory. Using the covariant matrix formalisms it is possible to calculate 
the entanglement at resonance of any trajectory analytically. Interestingly, through 
these resonances it is possible to produce motion generated quantum gates [19, 21, 29]. 
Recent results also show that genuine multipartite entanglement is generated through 
the periodic motion of the cavity [21]. 

The general case of a cavity undergoing an arbitrary acceleration and the possibility 
of performing actual experiments with mechanically oscillating optical cavities are under 
study [29]. In the cases of a linear sinusoidal or a uniform circular motion a resonance 
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appears at much lower frequencies for which no new photons are generated. The 
resonance leads to the generation of entanglement between existing and previously non- 
entangled cavity modes [29]. 

Understanding the effects of gravity on quantum properties and quantum 
information is of great interest. The analysis for entanglement in moving cavities already 
allows us to draw some interesting conclusions: via the equivalence principle [50] the 
results suggest that gravity can create entanglement in the field modes of a cavity (as 
long as the cavity is small). Consider a cavity in the vacuum held at a fixed position 
in a uniform gravitational field. The modes of the cavity would become entangled by 
changing the cavity's position and letting the cavity undergo periods of free fall. 

5.2. Localized wave packets 

An alternative for constructing localized field states is achieved by superimposing global 
plane wave modes Uk{x) = e^^'^^'' localized by a complete set of shape functions fm{k) via 
Umk{x) = J dkfjn{k)Uk{x) [41, 99]. A common example involves "box normalization" 
in which one quantizes in a box of length L and imposes periodic boundary conditions, 
leading to normalized states f^^ = /,-i/2g«27rma::/L_ ^^^^ ^^so normalize in both 

momentum and space degrees of freedom by using functions [41] fmi{k) = Xm{k)e~'^'^'^''^^'^ 
where Xm{k) = e~^/^ for \k — me\ < e/2. Here m, / are integers and e is a small positive 
constant with dimension of inverse length. The additional degree of freedom represented 
by the index / expresses the location of the box. In box normalization one takes the 
limit as L — 7- oo, while in the latter procedure one effectively uses and infinite number 
of boxes of length ~ 27re. 

5.3. Quantum communication with moving wave-packets 

As one recently proposed example of a relativistic quantum information protocol, 
Downes et al. [26] consider the case of Alice sending coherent states to Rob's accelerated 
detector, who then performs homodyne detection. Here the authors bypass the 
intermediate Unruh modes and consider the full integral Bogoliubov transformation (7) 
from which they form localized modes in frequency, as in the first equality in (61) in the 
appendix (section 8). In addition, they eschew the Schrodinger state transformation 
approach and instead consider the transformation of observables in the Heisenberg 
picture. With Alice stationary at an x-coordinate position to the right of focus of Rob's 
hyperbola at x = 1/a (x = 0), they consider only left moving L modes that travel along 
geodesies with x + t constant. Rob's homodyne operator is given by 0(r) = afj afj 6*"^ + 
h.c. where the superscripts K = S,L denote the signal and local oscillator modes (both 
of which are sent to Rob by Alice as localized Minkowski frequency modes) and the 
subscript i denotes modes (labeled by components of the wavevector), both parallel 
(longitudinal) and perpendicular (transverse) to Rob's acceleration. Rob integrates the 
photocurrent from his detector over a long time compared to the pulse length of Alice's 
signal. Thus the average value of the signal received by Rob is given by the expectation 
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value X = {f drOiir)) with variance V = {{J dTOi{T)Y) — X^, with the Minkowski 
coherent state = (a)L'f |0)m- Here Df{-i) = exp(7af - 7*af is the 

Minkowski displacement operator with 7 = a, /3 the amplitude of the signal and local 
oscillator respectively. The authors make reasonable approximations for matching the 
transverse spatial profile of the signal and local oscillator to Rob's detector ("beam- 
like" communication between Alice and Rob), and further assume that the longitudinal 
source wavevector is strongly peaked about a central frequency Ugo- They obtain 
the result that the normalized variance, defined as the variance V divided by the 
amplitude of the local oscillator as seen by Rob is given by the following Planck 
factor V = V/(5l = i/(e2'^^-(^+*) - 1). 

The "novel" dependence of the Planck factor \^ on a; + t is correctly attributed 
by the authors to the time dependent Doppler shift of the signal frequency Ugo sent 
by Alice as measured by Rob [42, 100, 101]. This should not be confused with the 
usual Planck factor l/(e^™'''« — 1) associated with the spectrum of frequencies w^.r 
Rob's detector measures when intercepting a Minkowski plane wave e'^^""^^'^^^ from 
Alice. The two are closely related as follows. The signal sent by Alice is centered 
on the 4- wavevector k^^ = (wso, — i^so, 0, 0), where we have taken k^^so = —<^so for 
a left moving photon. For Rob's trajectory = 1/a (sinh(ar), cosh(ar), 0, 0) his 4- 
velocity is given by = (cosh(ar), sinh(ar), 0, 0). By the previous discussion on 
Killing vectors in section 2.1.1, the frequency measured by Rob of this photon is given 
by ujr{t) = k^gU^^R = UsoG""^ ■ By substituting in x + t, the variance can be written as 
V = l/l^e^^'^^^'^^^'^ — 1), that is, the variance depends on the time dependent Doppler 
shifted frequency ujjiir) that Rob measures when he detects a Minkowski plane wave 
^iujso{x+t) _ gii^_R(r)/a q£ fixed Minkowski frequency Uso- This Minkowski plane wave is 
not detected as a positive frequency Rindler plane wave e"*^'*'^'^ by Rob, but rather as a 
combination of positive and negative frequency Rindler plane waves e^*'^'*-^'^. The square 
of the Fourier transform S of e*'^=°(^+*) as analyzed by Rob S = dr e""^'^-'^^ ^iujR{T)/a^ 
yields the Planck factor 1/ (e^'^'^''-^' — 1) [100, 101]. In fact, for finite integration 
times S{t) = JJ" (ir e"*'^'*'-^''^ e'"^^^''^/'* = JJ" dr e*'^^'^) the phase of the integrand has a 
stationary point [42, 102] at d(p{r)/dr = yielding the condition Ud^R = UsoG"''^ which 
is the frequency that appears in V. Previously wave packets for relativistic quantum 
information where considered in [103]. Recently a model on how to detect localized field 
modes has been introduced [27]. This model can also be used to study entanglement 
between localized wavepackets in non-inertal frames [28]. 

6. Change of state under Lorentz transformations 

While the main body of this work has been concerned with observer dependent 
entanglement for one or more observers undergoing constant acceleration, there is a 
large body of research on closely related effects for purely inertial observers (the zero 
acceleration case). This work is loosely classified under the heading of the change 
in quantum states under Lorentz transformations (LTs), and is intimately tied to the 
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concept of Wigner rotation and the unitary transformation of massive and massless 
particles under the Poincare group (translations, rotations and boosts). In the following, 
we review this research, which continues to be very active today, especially in connection 
to the transformation of entangled quantum states in curved spacetime. 

6.0.1. Preliminaries: Wigner rotation In flat (Minkowski) spacetime, the positive 
energy, single particle states of a massive particle forms a spinor representation of the 
inhomogeneous Lorentz (Poincare) group [104, 105] consisting of ten generators: four 
translations P^' = (P'^,P*), three generators {J*} of angular momentum that generate 
spatial rotations, and three generators {K^} of pure boosts. The quantum states, 
denoted by |p. A), are labeled by their spatial 3-momentum p (where = {p^,p) 
with p^ = E = ^/p"^ + ni? ) , and A the component of angular momentum along a 
quantization axis in its rest frame (typically taken to be along the third {z) spatial 
direction). Under a Lorentz transformation A the one-particle state transforms under 
the unitary transformation U (A) as [104, 105] 

f/(A)|p, A) = D'^AiW{^,p)) Ipa, A'), (44) 

A' 

where j is total angular momentum of the particle (equal to the spin of the particle 
in its rest frame), the summation is over A' = (— j, — j + l,...,j), and are the 
spatial components of the Lorentz transformed 4-momentum, i.e. p' where p'^ = h.^^p'^ . 
Most research has been primarily concerned with massive spin-| particles (j = ^), 
and massless spin-1 particles (j = 1, photons). While we will primarily discuss spin | 
particles, many related issues carry over to the case of photons, with important technical 
distinctions due to their absence of mass. In (44), D^l\(W{K,p)) is a (2j + 1) x (2j + 1) 
matrix spinor representation of the rotation group SU{2). The quantity W{K,p) is 
called the Wigner rotation (matrix), and depends on both the Lorentz transformation 
A describing the (constant velocity) inertial reference frame of the observer and the 
4-momentum p of the particle observed. The important point to note about (44) is that 
a single basis state |p. A) is transformed by a unitary representation of a momentum 
dependent rotation H^(A,p) into a superposition of all A (evaluated at the transformed 
momentum p^). 

The explicit form of the Wigner rotation matrix is given by 

W{^,p) = L"\Ap)AL{p), (45) 

where L{p) is a standard boost taking the standard rest frame 4-momentum k = 
(m, 0, 0, 0) to an arbitrary 4-momentum p, A is an arbitrary LT taking p — )■ Ap , and 
L~^{Ap) is an inverse standard boost taking the final 4-momentum Ap back to the 
particle's rest frame. Because of the form of the standard rest 4-momentum k, this final 
rest momentum k' can at most be a spatial rotation of the initial standard 4-momentum 
k, i.e. k' = W{A,p) k. The rotation group S0{3) is then said to form (Wigner's) little 
group for massive particles, i.e. the invariance group of the particle's rest 4-momentum. 



Observer dependent entanglement 



32 



The explicit form of the standard boost is given by [104] 

n(0 = ^ = coshe, 

U,{0 = ^ = smh^p\ L°=-sinhep., 

L'.iO = - (coshe - l)p'p„ 1,3 = (1,2,3), (46) 

where 7 = p^/m = E/m is the particle's energy per unit rest mass. Note that 
for the fiat spacetime metric rja/s =diag(l, — 1, — 1, — 1), po = p^ and pi = —p\ and 
P = P/\p\ = (sin 6' cos (f), sin 6 sin </>, cos 6) = is the direction of the 3-momentum p. 

In the rest frame of the particle, with (standard) 4-momentum k = (m, k = 0) the 
basis vectors |0, A) are eigenstates of the four momentum P^, total angular momentum 
= j (j + 1)) third component of angular momentum as (in units h = c = 1) 

P''\0,X) = P|0,A), 
J2|0,A) =j(j + l)|0,A), 

J^|0,A) =A|0,A). (47) 

The standard boost (46) can be decomposed [104, 105] as L{p) = R{p) Bs{\p\) R~^{p) 
where R{p) is a rotation that takes the z-axis to p, and i?3(|p|) is a pure boost in 
the z direction such that ||j?|i;,cr) = Ssdpl) |0, a) with |p| = m sinh^. The unitary 
representation of R{p) on the Hilbert space of states is given by U{R{p)) = e^'^'^'^ e*^"^^. 
Since the rightmost rotation in L{p) has no effect on |0, A) we have [105] 

|p. A) = L{p) |0, A) = R{p) BM) |0, A) = H{p) |0, A). (48) 

The states |p. A) are eigenstates of P^, J^, and the hehcity operator J ■ P/\P\ with 
eigenvalue A, the component of the particle's angular momentum along the direction p 
of it's linear momentum [105], 

P^|p,A) =p1p,A), 

J'IpA) =j{j + i)\p,X), 

^|p,A) = A|p,A). (49) 

(Note that J^\p,X) ^ A|p,A) since J now contains a contribution from the orbital 
motion). The Wigner rotation in (45) is given by W{A,p) = H~^{Ap) A H{p). The 
\p, A) are called helicity states. One can also define spin states denoted by \p, a) by 
suitable rotation |p, a) = D^^^{R^^{p))\p, A) [106], which are commonly used in the 
literature. Here R"^{p) rotates the 3-momentum direction p to the i:-axis. 

The formalism above carries through for massless particles, but with important 
physical differences. Since there is no rest frame for a massless particle, the standard 4- 
momentum is chosen to be /c = (1, 0, 0, 1) with k = and basis state |A;, A) labeled by .P 
as in (47) (which hold with substitution — )■ ). The general state |p. A) = H{p) A) 
is built up in exactly the same fashion as the massive case discussed in (48), where now 
the "energy" in i?3(|p|) is given by p^ = \p\ = e^. The helicity states \p, A) again satisfy 
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(49) with A now restricted to the two values ±j. The new feature for massless states is 
that the hehcity A is a Lorentz invariant since only the phase of the state \p, A) changes 
in a momentum dependent fashion under a LT A 

[/(A)|p,A) = e-^^^(^'P)|p,A). (50) 

Here, the notation indicates that 6{A,p) depends only upon the LT A and the direction 
p of the particle's momentum, but not upon its frequency The Wigner rotation angle 
0{A,p) can be computed as {k, X\W{A,p)\k, A), where again W{A,p) = H~^{Ap) A H{p) 
leaves the standard momentum k invariant kl = W{A^p) k (Wigner's little group is now 
IS0(2), rotations in plane perpendicular to k). For j = 1, 9{A,p) describes the amount 
of rotation experienced by the plane of polarization of a linearly polarized photon (for 
example calculations of 9{A,p) see [104, 107, 108]). 

Returning to the massive case, the spin states \p, a) (and similarly for the helicity 
states) are normalized as {p',a'\p,a) = {2tt)^ (2p^)5o-'aS^{p' ~ p)- One-particle states 
(wavepackets) are given by [105, 109] 



dM\Mp)\='^, (51) 



/oo PC 
dMMp)\p.<^). V^a(p) = / 

with respect to the Lorentz invariant measure dfi{p) = d^p/ {27r)^ i'^P^)- Under a LT A, 
the same state described in the boosted frame is given by 

/oo 
dMi;lip)\p,a), (52a) 

^Plip) = 5^Z}(f )(l^(A,A-V))V^<x(A-V) = (p,a|f/(A)|v[/), (525) 

V 

which can be obtained from a straight forward calculation using (51) and (44), the basis 
state normalization, a change of variables p — )■ A~^ p and utilizing dfi{A^^p) = dfi{p), 
along with some index relabeling. 



6. 0. 2. Early investigations Quantum information was first considered in a relativistic 
setting by Czachor in 1997 who analyzed a relativistic version of the Einstein- Podolsky- 
Rosen-Bohm experiment [107, 108, 109, 110, 111, 112, 113, 114]. Czachor's work showed 
that relativistic effects are relevant to the experiment where the degree of violation of 
Bell's inequalities depends on the velocity of the entangled particles. However, it was 
not until 2002 that the importance of investigating quantum information in the presence 
of spacetime was recognized. A. Peres & D. Terno pointed out that most concepts in 
quantum information theory may require a reassessment [109, 115]. Further interesting 
results on entanglement in flat spacetime were obtained by Alsing & Milburn [107], 
Adami, Bergou, & Gingrich [116, 117], and Solano & Pachos [118]. Their work showed 
that although entanglement is conserved overall under a change of inertial frame, it may 
swap between spin and position degrees of freedom. (For a comprehensive account of 
this early work up to 2004, with indications of future work on accelerated observers, see 
the review by A. Peres and D.R. Terno [109]). 
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The insightful and provoking work of Peres and Terno [109, 115] considered single 
particle wavepackets as in (51) with density matrix p = |^)(^| having components 
Paa'iP^P) — "^ctIp) V'cr' (p); ^^'^ traced out over the momentum to produce the reduced 
spin density matrix 

/oo 
dMrAp)i^Ap)- (53) 
-oo 

They showed that owing to the momentum dependent Wigner rotations that arise 
in the full state (52 a), that are induced by the LT A describing the reference frame 
of a new observer, the reduced quantity p^^"^ has no covariant transformation law, 
except in the limiting case of sharp momentum. Only the complete density p has a 
covariant transformation law. They concluded that the "spin state" of a single particle 
is meaningless if one does not specify its complete state, including momentum variables. 
Even though it may be possible to formally define spin in any Lorentz frame, there will 
be no relationship between the observable expectation values in different Lorentz frames 
[115]. 

Alsing and Milburn [107] and Terashima and Ueda [108] considered an EPR Bell 
state with factorable pure momentum eigenstates. In a second rest frame of another 
pair of observers A, B traveling perpendicularly to the initial rest frame, described 
by the pure boost A, the initial state |\E') is described as |\E'') = Uab{^) |^) = 
UA{^) ® Ub{A) I^). While the state |^ ) remains factorable between momentum and 
spin, the initial Bell state is coherently superposed with other Bell states. The reason 
is that Ua{-A) and Ub{A) give Wigner rotation angles of opposite signs for particles 
traveling in opposite directions. Maximal EPR correlations are recovered if the detectors 
in the boosted frame are rotated to the new rest frame direction of spin. 

Gingrich and Adami [116] considered wavepackets for two massive spin 1/2 particles, 
with similar form and transformation properties ([/ab(A) = Ua{-A)^Ub{A)) as (51) and 
(52 a), respectively 

I'^ab) = X] / J dMdM9arip,q)\p,(r)\q,T), (54a) 

err 
cr T 

Since LTs entangle the spin and momentum degrees of freedom within a single particle, 
entanglement can be transferred between them. Their work showed that this is also 
the case for pairs of particles, and that LTs affects the entanglement between the 
spins of different particles, as measured by the Wootters concurrence [119] on the 
reduced density matrix formed by tracing out over the momentum of the complete 
state pab = |^a_b)(^ab|- Thus, spin-spin entanglement is not a Lorentz invariant. 

These works generated an intense study of relativistic EPR correlations for both 
spin 1/2 particles(see [109, 120] and references therein) and photons [107, 109, 117, 121] 
that continues today. Much of the work for spin 1/2 particles (which we will primarily 
discuss henceforth) has focused on the relevant covariant observable (s) for spin such 
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that the expectation values obtained from measurements are the same in all inertial 
frames. Also in question is the meaning and validity of the reduced density matrix, 
especially in the case of tracing out the momentum from the complete quantum state. 
Terno et al. [115, 109] argued that in the relativistic setting momentum is a primary 
variable that has relativistic transformation laws that depends only on the LT matrix 
A that acts on the spacetime coordinates. Quantities such as spin and polarization are 
secondary variables that depend both on A and the primary variable momentum. Thus, 
even though the reduced density matrix for the secondary variable may be well defined in 
any coordinate system, it has no transformation law relating values in different Lorentz 
frames. Further, the unambiguous and seemingly natural construction of a reduced 
density matrix by means of tracing out over the primary variables is possible only if 
the secondary variables are unconstrained. In the absence of a general description, a 
case-by-case treatment is required. 

Shortly afterward, Czachor [122], Caban and Rembielinski [123] (see also [124] and 
references therein) showed that it is possible to define a Lorentz-covariant reduced spin 
density matrix for single massive particles. Such an object [123] contains information 
about the average polarization of the particle, as well as information about its average 
kinematical state. For sharp momentum, the reduced density matrix does not change 
under LTs. However, in the case of an arbitrary momentum distribution (51) the entropy 
of the reduced density matrix is in general not a Lorentz invariant. The conclusion of 
these studies is that while one can define a Lorentz-covariant finite dimensional matrix 
describing the polarization of a massive particle, in the relativistic case one cannot 
completely separate out kinematical degrees of freedom. These results ultimately stem 
from the momentum dependent spin transformations (44) induced by a LT. Extensions 
of this line of research have been carried out for states of two spin 1/2 particles (see 
e.g. [120] and references therein). Note that as pointed out by Peres and Terno [109] 
an invariant definition of entanglement of a pair of spin 1/2 particles typically utilizes 
the center-of-mass "rest frame" where (XI P) = 0- However, due to the problem of 
cluster decomposition [104] this definition is not adequate for more than two constituent 
particles since subsets of particles may have different rest frames. This is a difficult and 
still unresolved problem that has relevance to investigations of relativistic multipartite 
entanglement for more than two particles. 

6.0.3. Spin observable issues It is worth noting that complications in defining an 
appropriate reduced spin density matrix arise from the long appreciated non-unique 
definition of the "spin" operator for massive particles [125, 126, 127]. Ultimately, this is 
traced back to the fact that spin is only defined in the rest frame of the particle, where 
it coincides with the total angular momentum (since for p = the orbital angular 
momentum contribution vanishes). For studies of relativistic EPR correlations, most 
works have settled on some variation of the the spin operator defined from the Paul- 
Lubanski (PL) vector [127] W^j, = —l/2e^^ypaJ''^P'^ where J^^ are generators of the 
proper orthochronous Lorentz group with Jij = —Jji = ^ijkJk and Jjo = —Joi = Ki, 
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the generators of rotations and boosts, respectively. The two Casimir operators of the 
Poincare group are given by Ci = P^P^ = m"^ and C2 = W^W^ = —m'^j{j + 1) that 
commute with all the generators {P^, J^'^}, and label the basis states |p, cr) by their 
mass m and spin j. The PL vector = {J P, P^J + Px K) transforms as a 4- vector 
under LTs and is proportional to the total angular momentum in the rest frame of the 
particle W^f^g^ = (0, mJ). Hence, spin is well defined in the rest frame of the particle. In 
an arbitrary frame, the time component of the PL vector is proportional to the helicity 
operator J ■ P/\P\. 

For observers in an arbitrary inertial frame the spin S = J — L is defined as the 
difference between the total angular momentum, which is well defined as a generator 
of the Poincare group, and the orbital angular momentum L = R x P. While the 
momentum i-* is a well defined generator of the Poincare group, there is no well defined 
notion of a position operator (and hence the concept of localizability) in relativistic 
quantum mechanics. Hence L and S are not uniquely defined. Popular choices for the 
position operator [126] include the center-of-mass operator Rem. and Newton- Wigner 
position operator Rnw defined by 



7? -1 

-"-cm. 2 



P X K 

Rnw - Rem. - ^po(^ + po)- (55) 



This leads to the center-of-mass and Newton- Wigner spin operators 
„ W ^ 1 W^P 



^^w = — yW--— . (56) 

P^ m \ P^' + m / 

The operator Snw (favored by researchers such as Terno et al [115, 128] and Caban et 
al [123, 120]) satisfies the usual su{2) commutation relations and is the only axial-vector 
operator that is a linear combination of the PL vector. Using (55), a spin observable 
is defined by n ■ S]^w as the projection of spin along the unit spatial direction n. The 
downside of Snw is that by itself, it is not part of any known 4- vector, nor second 
rank tensor with well defined transformation properties under LTs. (Note that for both 
operators S ■ P = J ■ P so that we can also write = S ■ P). The center-of-mass spin 
operator Sc.m. (favored by researchers such as Czachor [110, 111, 115, 129, 128, 130]) 
does not generate the su{2) spin algebra and its eigenvalues are momentum dependent. 
Further, in contrast to the operator S%^r, the operator does not reduce to the 

second Casimir operator C2 = j(j + l). Despite these flaws, the spin observable n- Sc.m., 
has many favorable properties and is often used in the literature, as well as a normalized 
(nonlinear) version (introduced by Czachor [110]) Scz{'h) = fi ■ W /{m? + (n -p)^)^/^. 
A comparison of EPR correlations using both the Newton- Wigner and center-of-mass 
of spin observables [126] reveals that for fixed measurement directions, the EPR 
correlations do not necessarily decrease monotonically with velocity, and may exhibit 
local extrema for certain configurations (which also holds true for spin 1 particles). 
Furthermore, this effect occurs for both types of spin observables, and hence appears to 
be a generic feature of relativistic correlations, and in some cases strongly depends on 
which definition of spin observable is utilized. 
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Recent work by Friis et at [130] utilized Sc.m. for the investigation of a 
parameterized set of pure states of two spin 1/2 particles, with various configurations 
of entanglement between the spin and momentum degrees of freedom. As indicated by 
earlier works [107, 108, 110, 116] the maximum EPR correlation can be recovered in any 
inertial frame if both the state and the spin observable are Lorentz transformed. In an 
inertial frame in which the particle has 4-momentum p^^ = {p^,p), the authors show the 
observable defined by h{p) = {n- Sc.m.)/\^{n- Sc.m.)\, with X{n-Sc.m.) the eigenvector of 
the operator n ■ Sc.m., can also be written as fiiocaiip) = niocaiip) ■ Here cr is the usual 

3- vector of Pauli spin matrices, and riiocaiiji) = (\/l — /3^n_L + ^||)/(1 — /3^(1 — '^y))"'^^^ 
can be interpreted as the detector orientation riiocai as seen from the rest frame of the 
particle (where P = v/c is the velocity of the inertial frame). This later expression can 
also be written as hiocaiip) = [L~^{p)nf ai/\[L~^{p)ny\ where | [L~"^(p)n]-' | is the norm 
of the spatial portion of the Lorentz transformed orientation vector L~^{p)n. Here 
n = {n^, n) is the 4- vector whose spatial portion h is the orientation of the detector 
in the inertial frame with 4-momentum p. Thus, L~^{p) transforms the orientation 

4- vector n back to the rest frame uiocau whose spatial portion fiiocai indicates the 
local direction of the detector. Under a general LT A the orientation 4-vector is 
given by n" = An = AL{p)niocai so that n"(p) = [L^^(Ap)n"f ai/\[L~^(Ap)n"y\ = 
[W{A,p)niocaiy o'i/ \ [W {A, p)niocaiy \ - Since the Wigner rotation W{A,p) does not change 
the norm of hiocai we have that the observable h"{p) transforms under LTs as h"{p) = 
{W {A, p)niocai) ■ cr = U{A,p) {hiocai ■ U'^ {A,p). Thus, the expectation values of 
observable n"(j)) are the same as that in the rest frame of the particle niocai{p)- Most 
importantly, while the spin observable n"(p) depends on the momentum of the particle, 
the measurement direction n" corresponding to this observable does not [130]. 

Recently Saldanha and Vedral [124] have raised concerns about the physical 
implementations of the above spin observable measurements. They argue that while the 
above formulations are mathematically covariant, in the sense that spin measurements 
obtain values that are the same in all inertial frames, such measurements may be 
physically inconsistent. Consider a covariant description of the interaction Hpi of 
a measurement apparatus with spin observable S constructed from the PL vector 
W = {W°,W), where from (56) W^" = (5c.^.-P, P°5e.„i.) or W = (5wP,m5w + 
(P° — m){SNW ■ -P)-F*/I-Pp)- The interaction must be a Lorentz scalar of the form 
HpL = W^'G^ = W^G^ -W ■ G for some 4-vector = (G°, G) in order that 
the expectation values of a spin measurement are the same in all inertial reference 
frames. The authors state there are no known couplings of this type, and in particular, 
not for measurements where spin couples to the electromagnetic field. They argue, as 
have other researchers discussed above, that it is not possible to uniquely measure the 
spin of a particle independent from its momentum. They further conclude that any 
spin-momentum partition by means of reduced density matrices is "actually completely 
meaningless," while other researchers have argued that such bipartite partitions must 
be examined on a case by case basis, where meaningful information can be obtained. 
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6.0.4- Extension to curved spacetime The extension of the above work to curved 
spacetime for both spin 1/2 particles and photons has been examined by Terashima 
and Ueda [131], Alsing et al. [132, 133], Brodutch and Terno [134] and Palmer et al. 
[135]. In flat spacetime the inertial frame, described by the the LT A is global. In 
curved spacetime, this is true only locally at each spacetime point x = (x°,a;), where 
by the Equivalence Principle [38, 50], spacetime is locally (Lorentzian) fiat, and the 
rules of special relativity hold. The LT is now parameterized by the spacetime location 
X as local Lorentz transformation (LLT) which describes transformations 

between the instantaneous states of motion of different local observers at the same 
spacetime point x, for example freely falling, stationary (e.g. fixed radial coordinate), 
circular motion, or any trajectory with arbitrary acceleration. These LLTs are to 
be distinguished from general coordinate or world transformations dx'^/dx'^ which 
relates vectors, tensors, etc... describing the same spacetime quantity in different 
coordinates x and x'. The local Lorentz frame (LLP), i.e. the observer's local laboratory 
[136], has a small finite extent in space and time in as much as the curvature can 
be considered constant within this region. The LLP can be described by a tetrad 
e"„(x) [50, 137, 136], a set of 4 (orthonormal) local axes {e°(x)} labeled by (hatted 
local Lorentz) indices ((), 1,2,3) with spacetime components (world coordinate indices) 
a = (0, 1,2, ,3). The spatial triad {e*(x)} describes the orientation of the LLP in the 
surrounding spacetime, and the temporal axis e°(x), equal to the observer's (world) 
4- velocity (and thus tangent to their worldline), describes the observer's local proper 
time. As the observer moves through spacetime their LLP, described by their tetrad, 
twists and turns in the surrounding curved background. The observer describes events 
in their local laboratory by projecting world tensors onto this tetrad, e.g. a particle with 
momentum p°'{x) passing through the observer's local laboratory is described locally as 
p"(a:) = e"' ^{x) p^i^x). 

The quantum mechanical single particle states are now represented as \p'^{x),a) 
where p*(x) are the spatial components of the particle's 4- momentum p{x) = p"(x) ea(x) 
(also equal to p^{x)e^{x) where e^(x) can be taken as the coordinate basis vectors in 
the surrounding spacetime e^(a;) = (9^^). Terashima and Ueda [131] (see also appendix 
in Alsing et al. [132]) showed that as the observer with 4-momentum p{x) at x travels 
to a new, infinitesimally close spacetime point x' in proper time dr, the change 6p{x) = 
6p°'{x)ea{x) relative to the LLP can be written as a LLT 6p°'{x) = A"-(x) p''(x) (5r, 
so that A"-(a;) = + A"^(x) dr. One can now calculate the effect of the Wigner 
rotation expressed as a local version of (44) (with p — > p*(x), A — A(x), etc..) and 
(45) with an infinitesimal Wigner rotation given by W^^lx) = 6°"^ + i)"'^{x)dT, where 
the antisymmetric 'd°'^{x) only has non-zero components on the spatial indices (i.e. a 
rotation). Pinite Wigner rotations are obtained by a time-ordered integration along the 
observer's worldline with respect to his proper time. The implications of these results 
is that the orientation of the spin of the particle in the observer's LLP depends on 
instantaneous state of his motion as embodied in his tetrad e°'^{x) describing his local 
laboratory. This is not unexpected from the results of fiat spacetime investigation when 
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one invokes the Equivalence Principle. However, the entanglement of a pair of spin 
1/2 particles (initially created at a single spacetime point) when the parties are widely 
spacelike separated, depends not only upon the initial state of entanglement amongst the 
spin and momentum degrees of freedom (as in the fiat spacetime case), but also upon 
the full history of the motion through spacetime. That is, Wigner rotation induced 
spin-momentum entanglement occurs at each spacetime point along the trajectories of 
the constituent particles. While general relativity advocates an agnostic description 
of physics using arbitrary observers, the work of Alsing et al. [132] and Palmer et 
al. [135] shows that in the Fermi-Walker frame (the instantaneous, co-moving, non- 
rotating rest frame of the observer [50]) the Wigner rotation due to the effects of the 
gravitation field is zero. Though it appears contrary to the spirit of general relativity 
to single out a particular reference frame as special, the Fermi- Walker frame may prove 
advantageous to describe internal quantum mechanical degrees of freedom. While in 
the above discussion we have primarily concentrated on the Wigner rotation for massive 
particle in curved spacetime, many interesting results follow from the consideration of 
photons in an arbitrary gravitational field [133, 134, 135, 138, 139]. 

7. Open questions and future directions 

In prototypical quantum communication protocols Alice and Bob produce entangled 
resources (for example two entangled qubits) to be used to transit information. Alice 
and Bob take their systems to their labs. In case spacetime is fiat and Alice and Bob 
move at non-relativistic speeds the entanglement resources they produced will remain 
un-changed. However, what happens to quantum resources if Alice and Bob live in 
curved spacetime? If they move in the presence of a gravitational field? If they move 
at relativistic speeds? These are some of the most interesting questions in the field of 
relativistic quantum information. 

Understanding quantum information in relativistic settings might lead to finding 
ways to exploit relativity to improve quantum information tasks. A first step has been 
taken, finding suitable systems to store and process information in quantum field theory. 
The systems presented in this paper seam to be suitable candidates for this and hopefully 
ideas presented here will stimulate other researchers to find new systems and develop 
mathematical techniques that will allow us to find answers to the questions posed above. 
Quantum protocols such as teleportation are possible thanks to the tensor product 
structure of the Hilbert space. Constructing new protocols which are only possible in a 
quantum and relativistic world would indeed be exciting. 

An interesting question without doubt is if a covariant notion of entanglement 
can be constructed. The analysis of entanglement in relativistic settings has been 
performed by using the standard non-relativistic definition of entanglement and applying 
it to relativistic settings. By doing this we have found the entanglement is observer- 
dependent and that inconsistencies my arise when defining subsystems in special 
relativity. An open question is if there exits a more general notion of quantum 
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correlations in relativistic quantum theory which coincides with the standard notion 
in the non-relativistic regime. The answer to this question is intimately tied to issues 
related to the physical implementation of local measurements for observers in arbitrary 
motion. In this respect, the current research on a locally covariant formulation of 
quantum field theory in curved spacetime may prove useful [140, 141, 142]. 

Finally, we would like to say that we find the time ripe to find experimental 
demonstration of results in the field of relativistic quantum information. On one 
hand, cutting edge experiments in quantum information are approaching regimes where 
relativity starts playing a role and on the other hand, recent progress in the experimental 
demonstration of the dynamical Casimir effect [143, 144] might provide the requisite 
techniques to show that gravity has effects on entanglement . 
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8. Appendix 

A brief discussion of some of the states used in the literature to investigate 
acceleration/observer dependent entanglement is in order. A highly idealized scenario is 
to invoke the "single-mode approximation," (SMA) as considered by Alsing and Milburn 
[5], in which Minkowski annihilation operator is taken to be one of the right or left 
moving Unruh modes (since they do not couple to each other) in the integrand of (11) 

a-u^M ~ = cosh{rQ)anj - smh{rn)a'l jj. (57) 

One can then consider the entangled Minkowski state between Alice and Bob (left and 
right kets, respectively) 

1$) = i= (10^)-^ 10^,)-^ + lU^ 11.')-^) (58) 

when undergoes uniform acceleration (transforming into his Rindler alter-ego Rob) in 
wedge /. One can then expand Bob's Minkowski states |0(^/)-^ and ll^')"^ = al/ mI^'^')''^ 
in terms of the Minkowski vacuum expressed in terms of the Rindler wedge modes / 
and II given above and the adjoint of (57). This was done for idealized leaky cavities 
by Alsing and Milburn [5] and for free space modes by Fuentes and Mann [3] (see also 
the original Alsing and Milburn arxiv version in [5]). 
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Although the SMA gives degradation of entanglement with increasing acceleration, 
it is not rigourously justifiable since the Bogoliubov phase factors a^^^ ~ (w/)*'^^ = 
^itCiin{uii) g^j^g localized in frequency (the SMA was invoked to qualitatively represent 
the central frequency of some wavepacket). A more refined approximation was 
considered by Bruschi et al. [43] by considering the pre-accelerated initial Minkowski 
state to be composed of Unruh modes for Bob (compare with (58)) 

1$') = 4^ (|0^)^ lOf,)" + |ln)"). (59) 

V ^ 

The authors considered essentially a generalized version of the integrand in (11) by 
taking 

Oi^.M = qa an,R + qL^n^L, (60) 

with Igijp + Ig^^p = 1, which reduces to the SMA (57) for = 0. Again, entanglement 
was degraded as the acceleration is increased. 

The authors further considered the more physically realistic case of states with a 
finite spread of frequencies, i.e. wavepackets (see [41]), versus the highly idealized, 
sharply defined frequency states considered above (the latter of which capture the 
essence of the acceleration-dependent entanglement). Instead of (11) one considers 

/•CXD /»00 

«M = / dujf{uj)alM= dn[gR{n)a^^alj^ + gLin)a!^^alj= a^ji + a{, (61) 
Jo Jo 

where gR{fl) = du f {uj) a^^ and ghiS^) = du f {uj) a^^. By considering shape 

functions that are Gaussian in ln{ul) and adjusting parameters to ensure that the 

wavepacket version of the state (59) has negligible overlap between Alice's and Rob's 

state during acceleration, the SMA (57) is recovered with one of the qn, qi vanishing 

(see [43] for further details). 

The two global initial states 1$) (58) and |$') (59) consider different scenarios when 

the second party Bob undergoes uniform acceleration (to become the accelerated Rob). 

Both states |$) and |$') are maximally entangled, global Minkowski states. The first 

state 1$) is a "fixed" Minkowski state, independent of the acceleration a for both Alice 

and Bob, and is composed of the usual Minkowski plane wave u^j^m ~ g-^'^C'-'^^;) (^2). 

The second state |$') is composed of Minkowski states for Alice, and Unruh states for 

Bob whose mode u^^u is a linear combination of Rindler modes uq,j and that are 

both proportional to VL~^f'^{x — ety^^ in region / and II, respectively. On the one hand, 

by the transformation between Minkowski and Unruh modes (11), |$') Bob's portion of 

the state can be considered as an involved integration of global Minkowski modes. The 

entanglement in this state is then explored between the inertial Alice, and the Rindler 

observer Rob, who is specified by a particular value of the acceleration a, and whose 

mode function in region / is u^j. On the other hand, as discussed in section 3.1 on fiat 

spacetime entanglement, the state |$') may also be considered as a family of maximally 

entangled states parameterized by the dimensionless variable Q = u/a in non-inertial 

frames. By fixing the physical frequency u and changing a one analyzes the entanglement 

in a family of states which, all have the same frequency u as seen by observers with 
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different proper acceleration a. Alternatively, the entanglement corresponds to a family 
of states with different physical frequency u as seen by the same observer moving with 
fixed proper acceleration a [66]. 
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